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CHAPTER  1 


INTRODUCTION 

1.1  Problem  Definition 

National  interest  in  flight  vehicles  with  Mach  4  to  6  flight  capability  has  revitalized 
research  interest  in  high-speed  air  breathing  propulsion  systems  to  power  both  manned  and 
unmanned  vehicles.  With  increasing  flight  speeds  from  high  supersonic  to  low  hypersonic 
velocities,  the  entire  propulsion  system  will  operate  in  a  much  more  severe  environment.  As 
the  operational  limits  of  the  gas  turbine  engine  are  reached,  the  aerodynamic  conditions  will 
become  more  favorable  for  operation  of  a  subsonic-combustion  ramjet  propulsion  system. 
Thermal  management  of  propulsion  and  airframe  components  becomes  a  critical  design 
factor  beyond  the  Mach  6  flight  regime.  At  these  speeds,  freestream  air  cannot  be  used  for 
cooling  as  the  ram  air  temperature  will  approach  or  exceed  the  material  temperature  limits  of 
the  propulsion  system.  Therefore,  the  future  of  high-speed  flight  depends  to  a  large  extent  on 
the  ability  of  the  stored  fuel  to  be  used  as  a  heat  sink.  The  high  temperature  environment  also 
places  new  and  additional  constraints  on  the  combustion  processes.  For  example,  injection  of 
fuel  upstream  of  the  flameholder  is  no  longer  a  viable  design  option  for  mixing  and 
flameholding,  as  the  fuel  would  auto-ignite  in  the  high-temperature  inlet  air  and  destroy  the 
flameholder.  As  a  solution  to  this  problem,  an  alternative  to  a  flameholder  separate  from  the 
fuel  injector,  the  Integral  Fuel-Injector  Flameholder  (IFF),  was  studied  by  Hautman,  et  al. 
(1990).  The  IFF  also  serves  as  an  in-stream  fuel  injector  for  enhanced  mixing  and  could  be 
used  as  a  heat  exchanger  to  various  endothermic  type  fuels.  It  is  a  well-known  fact  that  fuel- 
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air  mixing  is  the  critical  factor  which  determines  the  combustion  efficiency  of  the  rambumer 
at  many  flight  conditions.  Previous  research  has  also  shown  that  the  use  of  vortex  generators 
attached  to  the  base  of  the  bluff  body  may  further  improve  mixing  and  combustion  efficiency 
much  like  lobed  mixers  used  with  turbofan  engines  (Spring  et  al.,  1992). 

The  flowfield  behind  geometries  such  as  bluff  bodies  has  been  the  topic  of  many 
studies  in  the  past.  Davies  and  Beer  (1971)  summarize  some  previous  experimental  work 
concerning  bluff  bodies.  It  is  also  well  established  that  the  recirculation  region  downstream 
of  such  bodies  enhances  fuel-air  mixing.  This  is  due  mainly  to  the  increase  in  local 
turbulence  levels  in  the  recirculation  region.  It  is  also  commonly  known  that  turbulence  is 
three-dimensional  as  well  as  rotational  (Tennekes  and  Lumley,  1972).  The  flowfield 
downstream  of  bluff  bodies  placed  in-stream  is  dominated  by  vortical  structures,  large  and 
small,  which  gives  the  flow  its  turbulent  characteristics. 

An  efficient  and  economical  combustor  design  methodology  requires  numerical 
experimentation  for  identifying  possible  prototype  candidates  before  fabrication  and  testing. 
Development  of  such  Computational  Fluid  Dynamics  (CFD)  code  requires  a  reliable 
database  against  which  the  predictions  can  be  compared  and  validated.  CFD  combustion 
engineers  will  require  a  detailed  database  on:  a)  the  combustor  entrance  conditions  such  as 
air  temperature,  pressure,  velocity,  and  turbulence  parameters  distributions  and  their 
relationship  to  inlet  operating  conditions;  b)  fuel  temperature,  composition,  and  distribution 
profiles;  and,  c)  location  of  the  flameholders,  and  injection  sites. 

The  flowfield  in  the  regions  immediately  downstream  of  such  bodies  being  three 
dimensional,  requires  a  three-dimensional  analysis.  The  introduction  of  laser  Doppler 
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velocimetry  (LDV)  has  allowed  for  the  nonintrasive  measurement  of  the  instantaneous 
velocities  and  other  turbulence  parameters  in  complex  turbulent  flows.  However,  the  use  of 
LDV  remained  limited  to  single  and  two-dimensional  studies  due  to  the  many  complications 
involved  in  measuring  the  third  component. 

1.2  The  Present  Research 

The  present  research  offers  a  detailed  study  of  the  flowfield  around  a  bluff  body 
mounted  in-stream  in  a  subsonic  isothermal  flow.  It  also  overcomes  the  difficulties  in 
obtaining  the  three  velocity  components.  A  benchmark  database  is  established  for 
development  and  validation  of  computational  codes.  In  order  to  establish  a  baseline  frame  of 
reference  for  future  work  employing  bluff  bodies  with  vortex  generators,  the  flow 
characteristics  of  a  cylindrical  bluff  body  are  experimentally  investigated  here.  Three- 
component  velocity  measurements  were  made  upstream,  and  in  the  highly  turbulent  flow 
downstream  of  the  bluff  body.  The  intense  wake  region  is  known  to  oscillate  at  a  frequency 
coupled  to  the  shedding  frequency  of  the  shear  layers  and  contains  large  negative  axial 
velocities.  Although  the  flow  behind  bluff  bodies  has  been  under  investigation  for  many 
years,  the  contribution  of  this  study  is  unique  as  all  three  velocity  components  were 
simultaneously  measured  at  each  point  without  introducing  disturbances  to  the  flow.  The 
three  velocity  components  were  measured  directly  using  three  orthogonal  and  overlapping 
probe  volumes.  This  gives  better  accuracy  than  the  more  common  method  of  obtaining  the 
third  velocity  component  as  described  by  Meyers  (1985)  and  Mehta  (1990).  This  more 
common  technique,  which  was  not  used  here,  uses  an  LDV  system,  which  allows  direct 
measurement  of  the  U  and  V  components,  together  with  a  second  LDV  system  rotated  at  a 
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certain  angle  to  allow  the  measurement  of  the  U  and  the  W  components.  Signals  from  the  U 
component  are  combined  with  signals  from  the  rotated  component  in  an  electronic  double 
balanced  mixer.  A  trigonometric  adjustment  of  the  difference  frequency  is  then  performed  to 
yield  the  W-component  measurement.  However,  if  the  W-velocity  component  is  small 
compared  to  the  other  velocity  components,  large  relative  errors  can  result. 

The  three  orthogonal  mean  velocity  components  (U,V,W)  were  measured, 
producing  profiles  in  the  transverse  and  spanwise  directions  at  1 1  axial  stations  in  the  far- 
field  and  at  8  axial  stations  in  the  near-fidld  behind  the  bluff  body.  Spanwise  velocity 
measurements  of  the  flowfield  immediately  downstream  of  the  bluff  body  were  also  taken. 
All  six  Reynolds  stresses,  normal  and  shear  (mm,  w,  ww,  uv,  uw,  vw),  and  all  nine  turbulent 
triple  products  (mmm,  mvv,  uww,  vmm,  wv,  ww,  wmm,  vwv,  www)  were  also  measured. 
Also,  this  report  includes  a  statistical  uncertainty  analysis  estimated  for  each  of  the  first, 
second  and  third  order  moments  using  Hrq  jackknife  method. 

Comparisons  between  commonly  used  approximations,  based  on  isotropic 
conditions,  for  vw,  mhw  and  vww  (when  the  third  component  is  not  available)  and  the 
experimental  measurements  are  made  and  presented.  Newly  derived  expressions  for  these 
parameters  are  suggested  using  the  direct  measurements  as  reference.  The  production, 
convection,  and  diffusion  of  turbulent  kinetic  energy  (TKE)  are  computed  directly  from  the 
experimental  data  using  numerical  differencing.  A  balance  of  TKE  in  the  flow  is  ^so 
performed  to  give  an  estimation  of  the  dissipation  of  TKE.  Regions  in  this  flowfield  where 
TKE  is  produced,  convected,  and  diffused  are  identified.  An  uncertainty  analysis  was  also 
estimated  for  each  of  the  TKE  terms.  Since  derivatives  with  respect  to  x  and  y  are  involved 
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in  the  TKE  balance  equation,  &st  order  derivatives  of  some  of  the  important  parameters  are 
also  reported  here.  Finally,  comparisons  between  some  recently  developed  numerical 
techniques  and  the  experimental  data  are  included. 

1.3  Background 

Simple  2-d  and  3-d  bluff  bodies  have  been  the  subject  of  experimental  investigations 
for  a  number  of  years  (Fuji  et  al.,  1978).  Bluff-body  flows  exist  in  many  places,  such  as  the 
flows  around  flameholders  in  gas  turbine  and  ramjet  engines,  and  flows  around  buildings, 
bridges,  and  support  structures.  Researchers  (Perry  and  Steiner,  1987;  Mansingh  and 
Oosthuizen,  1990)  have  found  that  the  turbulence  downstream  of  bluff  bodies  is  created  by  a 
regular  pattern  of  vortices  shed  from  the  trailing  edge  of  the  body,  and  that  symmetric  flow 
only  exists  as  a  time-average. 

1.3.1  Bluff  Body 

Some  examples  of  recent  bluff-body  studies  include  velocity  measurements  in  the 
flowfield  behind  flat  plates  by  Perry  and  Steiner  (1987).  They  used  a  flying-hot-wire  system 
for  all  measurements.  The  flat  plates  were  placed  normal  or  at  an  angle  of  attack  with  respect 
to  the  air  flow  at  Reynolds  numbers  that  ranged  from  14x10“  to  20x10“  based  on  the 
projected  height  of  the  plate.  Data  were  sampled  every  centimeter  (cm)  in  the  streamwise 
direction.  They  reported  that  the  flow  was  clearly  three-dimensional.  The  precise  nature  of 
the  three-dimensionality  was  not  stated.  However,  they  suspected  that  it  had  to  do  with  a 
periodic  waviness  of  the  vortex  cores  in  the  spanwise  direction.  They  also  concluded  that  the 
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mean-flow  patterns  of  the  bluff  plates  set  normal  to  the  freestream  flow  are  not  penGectly 
symmetrical,  and  are  highly  sensitive  to  the  end  conditions. 

Okamoto  (1981)  conducted  a  study  of  air  flow  past  a  low  aspect  ratio  hemispherical- 
cylinder  mounted  on  the  horizontal  surface,  in  a  wind  tunnel  with  a  400x400-inm  cross 
section  at  a  flow  Reynolds  number  of  4.74  x  10^  ,  based  on  the  diameter  of  the  cylinder.  By 
using  a  yaw  head  meter,  the  velocity  vector  was  obtained.  The  time-mean  velocity  and  static 
pressure  were  measured  by  pitot  and  static  pressure  tubes.  Also,  a  constant  temperature  hot 
wire  anemometer  was  used  to  measure  the  turbulence.  Distributions  in  the  turbulent  slrear 
flow  region  behind  the  hemispherical-cylinder  were  obtained.  These  measurements  showed 
that  the  recirculation  region  behind  the  bluff  body  extended  to  x/D  =  2.4.  This  recirculation 
zone  length  was  much  shorter  than  that  of  a  two-dimensional  fence  (x/D  =  13.2).  This 
difference  in  recirculation  zone  length  was  attributed  to  the  strong  downwash  behind  the 
hemisphere-cylinder.  A  Strouhal  number  of  0.152,  0.225,  and  0.222  was  reported  for  a 
hemisphere-cylinder,  a  circular  cylinder,  and  a  sphere,  respectively,  for  the  same  Reynolds 
number. 

Mansingh  and  Oosthuizen  (1990)  studied  the  flowfield  behind  a  rectangular  cross- 
section  cylinder.  They  determined  that  the  spanwise  base  pressure  measurements  at  mid¬ 
height  showed  considerable  three-dimensional  effects.  They  studied  the  effects  of 
introducing  a  splitter  plate  in  the  center  plane  of  the  rectangular  cylinder  at  various  dist^ces 
downstream.  Their  results  indicated  that  splitter  plates  altered  the  manner  of  vortex 
formation  in  the  wake  causing  a  decrease  in  shedding  frequency,  an  increase  in  base 
pressure,  and  a  reduction  in  the  overall  drag  by  up  to  50%.  Their  work  was  at  low  Reynolds 
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numbers,  typically  0.35x10^  and  1.15  xl0^  based  on  the  height  of  the  rectangular  cylinder 
and  the  freestream  velocity. 

Durao,  et  al.  (1988)  have  studied  the  flow  around  a  square  obstacle  mounted  in  a 
water  channel  at  a  Reynolds  number  of  1.4xl0\  based  on  the  obstacle  width  and  the 
upstream  velocity.  The  purpose  of  their  work  was  to  study  the  nature  of  the  fluctuations  in 
the  near-wake  region  of  two-dimensional  bodies.  Their  results  showed  that  in  the  zone  where 
the  shear  layer  oscillates  in  a  nearly  periodic  fashion,  the  turbulent  energy  associated  with 
the  random  turbulent  fluctuations  is  only  about  40%  of  the  total  turbulent  energy  calculated 
using  ensemble  averaging.  They  also  concluded  that  the  steady-state  numerical  calculation  of 
the  separated  flow  behind  two-dimensional  bodies  requires  modeling  of  the  periodic  vortex- 
shedding  motion. 

Williamson  (1989)  studied  the  three-dimensional  flowfield  in  the  wake  of  a  cylinder, 
mounted  vertically  in  a  water  tunnel.  Flow  visualization  was  carried  out  by  towing  a  cylinder 
through  water  in  an  xy  towing  tank,  using  fluorescent  dyes  washed  off  the  cylinder  surface, 
in  the  laminar  vortex  shedding  regime  (Reo<180),  a  discontinuity  in  the  Strouhal-Reynolds 
number  relationship  is  explained  by  a  transition  between  different  modes  of  oblique  vortex 
shedding.  Oblique  shedding  was  shown  to  be  caused  by  the  end  boundary  conditions  of  the 
cylinders  even  for  cylinders  hundreds  of  diameters  in  length.  They  found  that  by 
manipulating  the  end  conditions,  parallel  shedding  can  be  induced  which  results  in  a 
continuous  Strouhal  curve. 

Hautman  and  Rosfjord  (1989)  made  measurements  in  a  ramjet  test  combustor 
operating  at  realistic  flight  conditions.  Velocity  and  turbulence  intensity  measurements  were 
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made  using  a  two-channel  LDV  system  operated  in  the  back-scatter  mode.  High-speed 
photography  was  used  to  study  liquid  injection  and  atomization.  One  of  their  objectives  was 
to  assess  the  variation  of  inlet  turbulence  on  the  flowfield  downstream  of  a  bluff  body.  Two 
types  of  turbulence  modifiers  were  used  to  vary  the  upstream  turbulence.  Their 
measurements  showed  that  there  exists  regions  of  high  axial  turbulent  velocities  downstream 
of  the  Integral  Flameholder  Fuel  Injector  (IFF)  configurations.  Their  measurements  were 
then  used  for  the  development  of  ramjet  combustor  design  procedures. 

The  use  of  bodies  that  generate  vortical  and  turbulent  structures  has  been  suggested 
to  enhance  fuel  and  air  mixing  (Spring,  et  al.  1992).  It  is  believed  that  the  injected  fuel  from 
the  surface  of  the  in-stream  mounted  fuel  injector  will  follow  the  turbulent  flow  and  more 
thoroughly  mix  with  the  freestream  air.  McCormick  and  Bennett  (1993)  performed  a 
comprehensive  and  detailed  experimental  investigation  of  the  vortical  structure  in  the  free 
shear  layer  downstream  of  a  lobed  mixer.  They  showed  that  a  new  vortex  structure  exists  in 
addition  to  the  known  streamwise  vortex  array.  This  normal  vortex  sheds  periodically  from 
the  trailing  edge  and  plays  a  major  role  along  with  the  streamwise  vorticity  in  the  enhanced 
mixing  process. 

1.3.2  Laser  Doppler  Velodmetry 

In  recent  years,  laser  Doppler  velocimetry  (LDV)  has  become  increasingly  effective 
as  a  tool  for  non-intrusively  measuring  the  instantaneous  velocities  at  points  in  complex 
turbulent  flows.  This  is  because  LDV  is  a  reliable  instrument  for  making  instantaneous 
velocity  measurements  in  highly  turbulent  flows  and  in  flows  with  separation  and 
recirculation  regions.  However,  few  simultaneous  three-component  measurements  have  been 
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reported  to-date,  due  mainly  to  the  increase  in  instrumentation  and  signal  processing 
equipment  required  to  make  these  measurements.  In  addition,  optical  access  is  required  for 
all  three  components.  Other  difficulties  include  the  ability  to  maintain  a  high  coincident  data 
validation  rate  so  that  temporal  flow  characteristics  can  be  resolved.  Also,  the  continuous 
drift  (due  to  vibration  and  room  air  temperature  fluctuations,  etc.)  of  the  three  individual 
probe  volumes  necessitates  adjustment  for  maintaining  coincident  three-component  data. 

Sjunnesson,  et  al.  (1991)  used  a  two-component  LDV  system  to  study  the  flowfield 
behind  2-d  bluff  bodies  of  various  shapes  at  Reynolds  numbers,  based  on  the  flameholder 
height,  ranging  from  2.8x10“  to  8.4x10“  for  the  isothermal  and  reacting  flow  conditions  in 
air.  Their  study  revealed  that  the  recirculation  zone  was  longer  for  the  reacting  flow  case, 
that  the  flow  recovered  in  a  shorter  distance  in  the  isothermal  case,  and  that  the  strong  vortex 
shedding  present  in  isothermal  flows  was  not  found  in  the  reacting  flow. 

Geropp  and  Leder  (1985)  used  a  single  component  (LDV)  to  document  the  flowfield 
behind  a  2-d  wing  body.  They  changed  the  angle  of  attack  of  an  airfoil  and  a  wedge,  and 
showed  that  the  length  of  the  recirculating  region  in  the  wake  increases  for  all  bodies  with 
increasing  angle  of  attack.  The  suppression  of  vortex  interaction  was  also  obtained  by 
placing  a  splitter  plate  at  mid-height,  downstream  of  the  body,  which  resulted  in  a  reduction 
in  drag. 

The  turbulent  wake  behind  a  circular  cylinder  was  .studied  by  Absil  (1988)  u.sing  a 
two-color  dual-beam  laser  Doppler  anemometer  (LDA).  The  Reynolds  number  based  on  the 
diameter  of  the  cylinder  was  1300.  Mean  velocities,  turbulence  intensities,  Reynolds  stresses, 
skewness  and  kurtosis  factors,  and  other  correlation  functions  were  reported.  His 
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measurements  consisted  mainly  of  stations  in  the  far  wake  where  no  reverse  flow  existed. 
This  was  done  on  purpose  to  avoid  the  use  of  a  Bragg  cell  in  his  LDA  system.  His  data 
showed  reasonable  agreement  with  hot  wire  data  taken  previously.  The  investigation 
concluded  by  demonstrating  that  LDA  is  a  powerful  tool  for  measuring  various  turbulence 
quantities  in  turbulent  flows. 

Hosokawa,  et  al.  (1993)  studied  the  effects  of  flameholder  geometries  on  the 
flowfield  characteristics  such  as  the  recirculation  zone  size  and  vortex  scale.  They  used  a  Ld 
fiber-optic  based  laser  Doppler  velocimeter  (FLDV)  in  an  isothermal  air  flow  behind  two 
different  versions  of  V-gutter  shaped  flameholders,  with  and  without  an  extended  base.  They 
measured  the  axial  mean  and  fluctuating  Velocities,  and  the  frequency  of  the  fluctuating 
flow.  They  found  that  the  rectangular  section  of  the  extended  flameholder  played  a  role  in 
the  flow  stratification,  due  to  the  uniform  velocity  profile  over  the  extended  surface.  This 
also  increased  the  fluctuating  frequency  and  decreased  the  length  of  the  recirculation  zone. 

Larousse,  et  al.  (1991)  used  a  two-dimensional  LDV  system  to  study  the  turbulent 
flow  behavior  behind  2-d  and  3-d  surface-mounted  obstacles  in  a  channel,  where  air  flowed 
at  a  Reynolds  number  of  lO’,  based  on  the  obstacle  height.  Major  differences  in  the  flow 
patterns,  both  upstream  and  downstream  of  the  obstacles,  were  noted.  A  distinct  unsteadiness 
of  the  upstream  flowfield  for  three-dimensional  obstacles  led  to  increased  values  in  the 
Reynolds  shear  stress. 

1.3.3  Computational  Fluid  Dynamics 

Durst  and  Rastogi  (1977)  used  a  one-component  LDV  system  to  study  the  flowfield 
behind  a  square  obstacle  in  a  two-dimensional  water  channel.  The  measurements  were  made 
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at  a  Reynolds  number  of  1.7xl0^  Theoretical  investigations  were  also  made  using  finite 
difference  calculation  procedures  for  recirculating  flows  and  boundary  layers  in  order  to 
solve  the  governing  fluid  flow  equations.  The  turbulent  viscosity  was  obtained  from  the 
solution  of  the  two-equation  k-8  turbulence  model.  For  boundary  layer  regions,  the  three- 
equation  k-E-Mv  turbulence  model  was  also  used.  They  found  that  the  calculated  longitudinal 
velocities  did  not  compare  well  with  those  measured.  They  suggested  that  the  k-£  turbulence 
model  needs  to  be  modified  in  order  to  obtain  more  accurate  calculations  in  the  separated 
flow  regions. 

Castro  (1977)  studied  the  numerical  accuracy  of  finite  difference  techniques  typically 
used  for  the  calculation  of  elliptic  type  turbulent  flows.  The  size  of  the  truncation  error  term 
was  investigated  both  analytically  and  by  direct  comparison  between  predictions  and 
experiments.  First,  a  classical  turbulent  shear  flow  was  considered.  For  this  flow,  it  was 
found  that  although  the  local  numerical  errors  became  large,  the  general  features  of  the  flow 
were  well  predicted.  Next  the  flow  over  a  backward-facing  step  and  over  a  two-dimensional 
fence  were  considered.  It  was  concluded  that  inadequate  turbulence  modeling  was  not  the 
sole,  or  at  least,  not  the  major  cause  of  discrepancies  between  predictions  and  experiments.  It 
was  further  observed  that  common  numerical  methods  with  typical  grid  densities  are  often 
incapable  of  accurate  predictions  due  to  severe  numerical  errors  in  the  regions  of  difficult 
geometry  such  as  sharp  corners. 

Moss,  et  al.  (1977)  used  pulsed-wire  anemometry  to  measure  mean  velocities  and 
Reynolds  stresses  in  turbulent  (Re  =  SxlO'')  recirculating  flow  past  a  number  of  bluff  bodies. 
The  distributions  of  the  Reynolds  stresses  were  given  for  the  recirculating  zones  and  the 
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mixing  layers.  Measurements  of  the  mean  velocity  and  turbulence  energy  were  compared 
with  the  predictions  of  a  computer  program  based  on  a  two-equation  turbulence  model.  - 

Przulj  and  Younis  (1993)  performed  a  numerical  simulation  study  of  turbulent  vortex 
shedding  for  flow  over  a  square  cylinder,  at  a  Reynolds  number  of  250  based  on  the  cylinder 
height.  Their  model  captured  the  experimentally  documented  size  of  the  recirculation  zone. 
They  showed  that  the  prediction  of  turbulent  vortex  shedding  is  possible  with  a  standard  k-e 
turbulence  model.  However,  the  model  was  found  to  lack  sensitivity  to  the  imposed  mean- 
flow  unsteadiness.  They  presented  an  assessment  of  important  factors  such  as;  a) 
discretization  of  convection  terms;  b)  the  use  of  the  log-law  for  simulation  of  the  near-wall 
regions;  and  c)  the  choice  of  turbulence  model  which  influences  the  time-averaged 
simulations  of  turbulent  vortex-shedding  associated  with  bluff  bodies.  They  suggested  that 
suppression  of  vortex-shedding,  which  has  often  been  associated  with  the  use  of  k-£ 
turbulence  modeling  was  due  to  inappropriate  numerical  discretization  schemes.  Turbulence 
models  that  have  been  developed  for  steady  flow  are  shown  to  be  inadequate  for  vortex- 
shedding  applications.  A  simple  modification  was  proposed  for  sensitizing  the  standard 
turbulence  models  for  such  flows.  After  modification  of  the  £-production  term,  the  model 
accurately  predicted  the  measured  axial  mean-flow  and  turbulence  parameters.  No  transverse 
mean-flow  velocities,  or  turbulent  stresses  were  reported. 

Martensson,  et  al.  (1991)  used  different  numerical  simulation  techniques  to  predict 
the  flow  behavior  downstream  of  two  different  bluff-body  geometries,  mounted  in-stream. 
The  Reynolds  numbers  investigated  were  4.5x10“  and  8.7xl0“,  based  on  the  bluff-body 
height.  A  triangular  cross  section  geometry  and  a  triangular  cross  section  geometry  with  a 
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square  cross  section  afterbody  were  studied.  The  bluff  bodies  were  mounted  with  the  apex  of 
the  triangular  section  facing  the  air  flow.  They  reported  results  from  a  number  of  numerical 
experiments,  using  different  grids,  numerical  schemes  and  flow  conditions.  The  3-d 
numerical  results  generally  agreed  better  with  experimental  results  than  did  the  2-d 
numerical  results.  There  was  a  major  improvement  in  matching  the  transverse  turbulent 
fluctuating  term,  v',  as  well  as  the  Reynolds  shear  stress,  UV,  with  the  experimental 
measurements,  when  using  3-d  code.  The  spanwise  velocities  were  not  measured. 

Kodali  and  Amano  (1986)  conducted  a  numerical  study  on  the  turbulent  flowfield 
behind  bluff  bodies.  They  reported  on  an  external  flow  over  a  disk  and  an  internal  flow  past  a 
cone.  They  compared  two  k-e  turbulence  models,  and  two  Algebraic  Stress  Models  (ASM) 
where  the  only  differences  were  the  values  of  modeling  constants.  When  compared  with 
actual  measurements,  their  results  indicated  that  the  ASM  generally  gave  better  predictions 
in  the  recirculating  region  than  did  the  k-£  model.  They  also  found  that  a  modification  to  the 
ASM  model  gave  a  significant  improvement  of  prediction. 

Franke  and  Rodi  (1991)  calculated  the  vortex  shedding  behind  a  surface  mounted 
square  cylinder  at  a  Reynolds  number  of  2.2x1 0^  Three  variants  of  the  k-e  eddy-viscosity 
model  and  two  Reynolds-stress  equation  models  were  used  to  model  the  turbulence 
characteristics  of  the  flow.  The  experimental  data  used  were  from  Lyn  (1989)  who  used  a 
two-component  LDV  system,  and  from  Durao,  et  al.  (1988)  for  a  Reynolds  number  of 
1.4xl0\  They  found  that  the  variations  between  the  k-e  models  and  the  Reynolds  Stress 
models  (RSM)  were  due  to  the  near-wall  treatments.  They  concluded  that  the  RSM  models 
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do  a  better  job  in  predicting  the  length  of  the  recirculation  region  and  the  level  of'  total 
turbulence  fluctuations. 

More  recently,  numerical  studies  using  various  CFD  techniques  continued  to  be 
applied  to  simple  bluff-body  geometries.  These  studies  have  produced  solutions  using 
techniques  varying  in  sophistication  from  steady-state  algebraic  stress  models,  to  the 
dynamic  subgrid-scale,  Large-Eddy  Simulations  (LES)  of  Yang  and  Ferziger  ( 1993). 

Some  early  solutions  modeled  only  half  of  the  geometry,  thus  imposing  synmietry  on 
the  flow.  Time-averaged  measurements  indicated  that  the  flowfield  was  symmetric  with  a 
steady  pair  of  vortices  downstream  of  the  bluff  body.  However,  detailed  studies  have  shown 
that  the  instantaneous  flowfield  consists  of  a  structured  pattern  of  alternating  large  scale 
vortices.  A  symmetry  plane  enforced  by  the  numerical  model  precludes  the  formation  of-the 
shed  vortices,  and  thus  requires  the  turbulence  model  to  predict  the  proper  turbulent 
fluctuations  in  the  wake.  Comparisons  to  experimental  measurements  have  shown  such 
predictions  to  be  inadequate  for  bluff-body  flowfields. 

The  LES  technique  is  formulated  and  applied  to  a  plane  channel  flow  in  which  a  two- 
dimensional  obstacle  is  mounted  to  the  tunnel-surface  as  reported  by  Yang  and  Ferziger 
(1993).  The  Reynolds  number  was  3.2xl0\  based  on  the  height  of  the  obstacle.  A  Direct 
Numerical  Simulation  (DNS)  was  also  conducted  for  comparison.  They  noticed  that 
although  the  obstacle  is  two  dimensional,  the  instantaneous  separation  and  reattachment  lines 
were  found  to  be  far  from  two  dimensional.  Secondary  and  tertiary  unsteady  flow  regions 
were  found  near  the  bluff  body.  No  comparisons  with  experimental  measurements  were 
made. 
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Wemer  and  Wengle  (1991)  applied  the  LES  technique  to  model  the  flow  downstream 
of  a  3-d  bluff  body  (cube)  mounted  on  the  surface  of  a  plate  channel,  at  a  Reynolds  number 
of  5x1 0\  After  comparison  of  their  predicted  results  with  the  experimental  data  of  Larousse, 
et  al.  (1991),  they  recommended  that  an  improved  spatial  resolution  immediately  above  the 
top  face  of  the  cube  and  near  the  two  sides  was  necessary.  In  the  case  of  the  shear  stress 
comparison,  they  concluded  that  second-order  statistical  moments  (like  mv)  were  not 
sufficiently  smooth  in  certain  regions  of  the  flow,  such  as  just  downstream  of  the  bluff  body, 
and  at  the  reattachment  region.  This  was  attributed  to  the  slowly  converging  second-order 
statistics  in  these  particular  flow  regions. 

The  Renormalized  Group  (RNG)  theory  was  applied  by  Choudhury,  et  al.  (1993)  to  a 
selected  number  of  separated  turbulent  flows,  such  as  a  backward-facing  step  and  the 
periodic  vortex  shedding  behind  a  circular  cylinder.  A  substantial  improvement  over 
conventional  turbulence  models  was  realized,  especially  in  predicting  the  mean  flow 
quantities  and  in  capturing  the  vortex  shedding  frequency,  when  the  RNG  k-e  model  was 
used. 

Also,  Yakhot  and  Orszag  (1992)  reported  excellent  results  for  homogeneous  shear 
flow  and  for  flow  over  a  backward-facing  step.  For  the  shear  flow,  the  modified  Reynolds 
stress  transport  model  was  also  used.  The  modified  eddy  viscosity  models  gave  the  best 
results  for  the  backward-facing  step  flow.  In  all  cases  no  adjustments  to  the  model  turbulence 
constants  were  made. 
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1.4  Outline 


The  primary  objective  of  this  research  was  to  perform  an  extensive  experimental 
study  of  the  flowfield  downstream  of  a  2-d  bluff  body,  mounted  at  mid-stream  iit  a 
rectangular  cross-section  test  section.  The  test  conditions  are  atmospheric  pressure  and 
temperature.  An  experimental  3-d  mapping  of  the  flowfield  in  the  near-  and  far-field  regions 
downstream  of  as  well  as  upstream  of  the  bluff  body  is  presented.  The  velocities  were 
measured  with  a  state-of-the-art  3-d  LDV  system.  This  allowed  for  nonintrusive  pointwise 
measurements  throughout  the  tunnel  test  section.  The  flowfield  properties  of  interest  include 
mean  and  fluctuating  axial,  transverse,  and  spanwise  velocities,  turbulent  normal  and 
Reynolds  stresses  and  turbulent  triple  products.  Comparisons  of  some  of  these  measurements 
with  CFD  generated  predictions  are  also  presented. 

Chapter  2  is  a  review  of  turbulence  modeling  and  is  presented  primarily  to  show 
which  turbulence  parameters  are  of  interest  to  turbulence  modelers.  This  section  is  offered  as 
guidance  to  the  experimentalist  so  that  useful  turbulence  measurements  result. 

The  experimental  apparatus,  designed  and  used  for  this  study,  is  described  in  detail  in 
Chapter  3.  The  LDV  system,  flow  system,  the  bluff-body  geometry,  and  the  test  section 
design  are  presented  in  this  chapter. 

Chapter  4  describes  the  test  procedure  and  summarizes  the  test  run  conditions.  A 
discussion  of  measurement  uncertainty  and  errors  is  also  given  here. 

The  far-field  velocity  measurements  and  turbulence  parameters  are  presented  in 
Chapter  5.  The  axial  stations  include  two  upstream  and  nine  downstream  locations. 
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The  near-field  velocity  measurements  and  turbulence  parameters  are  presented  in 
Chapter  6.  Also  included  are  experimentally  approximated  expressions  for  and 

uww.  A  turbulent  kinetic  energy  balance,  used  to  estimate  the  turbulent  dissipation  rate,  is 
also  presented  using  the  measurements  made  in  the  near-field  region.  A  statistical 
uncertainty  estimate  for  each  of  the  TKE  terms,  at  the  recirculating  bubble  station  location 
(x/H=0.9),  is  also  presented.  Spatial  gradients  of  some  important  fluid  parameters  and  a 
check  of  the  two-dimensionality  of  the  flow  are  also  included. 

Comparisons  of  the  numerical  predictions  with  experimental  measurements  are 
presented  in  Chapter  7. 

Finally,  conclusions  and  recommendations  are  given  in  Chapter  8. 


CHAPTER  2 


TURBULENCE  MODELING 

2.1  Introduction 

Most  flows  occurring  in  the  area  of  fluid  mechanics  and  related  applications  are 
turbulent  in  nature.  No  specific  definition  can  be  given  to  such  a  phenomenon.  However, 
turbulence  is  commonly  described  by  its  flow  characteristics  (Tennekes  and  Lumley,  1972). 
Turbulent  flows  are  random  and  have  unique  characteristics;  every  turbulent  flow  is 
different.  In  general,  turbulent  flows  exhibit  large  irregular  spatial  and  temporal  fluctuations 
in  velocity,  pressure,  temperature,  and  vorticity.  It  is  this  irregularity  of  the  flowfield  that 
necessitates  the  usage  of  statistical  averaging  methods  for  the  determination  of  turbulence 
parameters,  a  technique  widely  accepted  in  engineering.  Because  turbulent  flows  are  so 
prevalent  and  interdisciplinary,  there  has  been  an  increased  interest  in  better  understanding 
turbulent  flow  phenomena.  In  particular.  Computational  Fluid  Dynamics  (CFD)  has  surged 
as  a  new  and  very  important  methodology  for  predicting  and  better  understanding  the 
turbulent  flow  phenomena.  One  driving  force  for  this  activity  in  turbulence  modeling  is  the 
recent  development  of  high-speed  digital  computers,  needed  to  solve  the  mathematical 
equations  that  attempt  to  model  complex  problems  such  as  those  encountered  in  turbulent 
flows.  One  of  the  underlying  reasons  for  the  interest  of  making  the  measurements  reported  in 
this  study,  is  the  validation  of  the  computational  approach  (or  numerical  code)  for  predicting 
the  mixing  behind  bluff-body  flameholders.  So  far,  most  practical  turbulence  models  fail  to 
properly  predict  or  capture  the  basic  features  of  turbulent  separated  flows  such  as  those 
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resulting  from  vortex  shedding  behind  a  bluff  body.  The  nature  of  complex  flows  demand 
the  solution  to  the  nonlinear,  rotational,  three-dimensional  and  time-dependent  flow 
equations.  Some  of  the  parameters  of  interest  are  the  magnitude  of  the  velocity  vector, 
turbulence  intensities,  the  Reynolds  stresses,  the  vortex-shedding  frequency  and  the  length  of 
the  recirculation  region.  A  brief  description  of  turbulence  modeling  is  presented  next. 


2.2  Ori2in  and  Theory  of  Turbulence  Modeling 

The  subject  of  turbulence  began  to  be  better  understood  some  100  years  ago,  when 
Osborne  Reynolds  (1895)  came-up  with  the  idea  of  time-averaging  the  Navier-Stokes 
equations  (NSE).  Since  turbulence  is  a  continuum  phenomenon,  the  time-dependent,  three- 
dimensional  NSE  should  give  a  sufficient  description  of  the  physics  of  turbulent  flow. 
Consequently,  statistical  averaging  methods  using  time-averaging  must  be  used,  which  in  turn 
lead  to  statistical  correlations  which  must  be  determined  in  order  for  closure  to  take  place.  The 
dynamic  behavior  of  laminar,  as  well  as  turbulent  flows,  for  a  Newtonian  fluid  in  the  absence 
of  any  external  body  forces  is  governed  by  the  momentum  transport  balance,  written  in 


Cartesian-tensor  notation  as; 

a(pt/,) .  a(pf/,f/,)  dp  a  ^  ac/,  ^ 

dt  Bxj  ax,/a;c/^^ajc/ax, 


(2.1) 


The  conservation  of  mass  principle  necessitates  the  use  of  the  continuity  equation  as; 

£;(pt',)  =  0  (2.2) 


Equation  (2.2)  together  with  the  momentum  transport  balance  (2.1)  constitute  the 


Navier-Stokes  equations.  These  equations  are  given  in  rectangular  coordinates  here  so  as  to  be 
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useful  for  the  eurrent  study.  The  instantaneous  velocities  U.  are  then  decomposed  into  a  mean 
velocity  U  ■, ,  and  velocity  fluctuations  m,  ,  such  that; 

U.=Ui+u.  (2.3-) 


Where  f/,  ,  is  interpreted  as  the  time-averaged  velocity,  defined  as; 


Ui=  —  ]U,dt  (asT->oo) 


(2.4) 


-T 


(The  mean  value  of  a  fluctuating  quantity  is  zero  by  definition.) 

These  definitions  are  then  substituted  into  the  Navier-Stokes  equations,  and  expanded 
giving  the  time-averaged  momentum  equations  for  turbulent  flow  also  called  the  Reynolds 
averaged  equations.  The  governing  equation  is  given  in  tensor  form  as; 


d(pui)  ^dipUiUj)  ap.  a  .  .at/,  ,  at/^  — , 

- + ^ +  +  3— )- 


at 


dxj 


Bxi  dXj  dXj  dx, 


(2.5) 


This  equation  is  identical  to  the  NSE  (2.1)  except  for  the  appearance  of  second-moment 
correlation  involving  the  turbulent  velocity  fluctuation  terms,  Similar  to  the  pressure 

and  viscous  stresses,  these  fluctuating  inertia  terms  act  on  an  element  of  the  fluid  and  have  the 
dimensions  of  shear  stress  (force/area).  These  are  turbulent  stresses  induced  into  the  flow 
because  of  turbulence.  The  additional  stresses  are  identified  after  their  originator,  Reynolds, 
and  are  commonly  known  as  the  “apparent”  or  “turbulent”  or  “Reynolds”  stresses.  They  add 
complexity  to  the  problem  by  introducing  six  additional  equations  to  the  three-dimensional 
flow  equations,  which  must  be  prescribed  first  in  order  to  solve  Equation  (2.5)  for  the  mean- 
flow  parameters.  Since  we  have  more  unknowns  than  we  have  equations,  the  objective  of 
turbulence  modeling  is  to  define  additional  equations  which  relate  these  unknown  Reynolds 
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stresses  to  the  mean-flow  parameters.  Since,  these  equations  cannot  be  determined  a  priori,  this 
process  is  known  as  turbulence  closure. 

Turbulence  modeling  was  initially  formed  on  the  mathematical  derivation  of  the 
turbulent  stresses  using  the  same  procedure  used  for  laminar  flows.  The  predominant 
assumption  was  that  the  effects  of  turbulence  in  the  flowfield  can  be  related  to  mean-flow 
parameters,  such  as  the  mean  velocity  gradient.  Boussinesq  (1877)  was  the  first  to  propose  the 


concept  of  an  “apparent”  or  “turbulent”  or  “eddy”  viscosity.  In  order  to  account  for  the 
presence  of  the  Reynolds  stresses  in  the  turbulent  momentum  balance  of  Equation  (2.5),  it 
becomes  necessary  to  model  and  evaluate  these  Reynolds  stresses  suitably.  Boussinesq 
reasoned  that  the  apparent  turbulent  shearing  stresses  act  similar  to  viscous  stresses  and  thus 
might  be  related  to  the  local  rate  of  mean  strain  through  an  apparent  turbulent  viscosity  (see 
Hinze,  1975).  He  formulated  the  concept  of  an  eddy-viscosity  defined  as  follows; 


-  dUi  dUj 


(2.6) 


The  proportionality  factor,  p,,  is  known  as  the  “apparent  turbulent  coefficient  of 
viscosity”  or  “eddy-viscosity”  and  is  analogous  to  the  constant  of  proportionality  of  the  fluid’s 
laminar  molecular  viscosity,  p.  The  problem  now  is  to  find  a  suitable  value  or  relation  for 
determining  p,.  Boussinesq  originally  assumed  a  constant  value  of  eddy-viscosity,  p,. 
However,  unlike  laminar  flows,  turbulent  flows  are  nonhomogeneous,  and  a  single  value  for  p, 
would  imply  homogeneity.  The  ‘turbulent  viscosity’  of  the  fluid,  p„  unlike  the  laminar 
viscosity  which  is  a  property  of  the  fluid,  is  a  flow-dependent  variable  since  turbulence  is  a 
characteristic  of  the  flow.  Another  problem  lies  in  the  mathematical  contraction  of  Equation 
(2.6).  If  i  is  set  toj  and  summed  from  i  =  i  to  3,  the  right-hand  side  is  zero  as  required  by 
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continuity,  but  the  left-hand  side  is  zero  only  if  there  is  no  turbulence.  Therefore,  to  correct  for 
this  deficiency.  Equation  (2.6)  is  modified  by  an  additional  term  as  shown  by; 


-  ^dUi  2^  , 

(2.7) 

Where  k  is  the  turbulent  kinetic  energy  (TKE),  and  5,^  is  the  Kronecker  delta  (5,^.  =  7  if  i  =j 
and  d y=  0  if  i^f). 

Substitution  of  Equation  (2.7)  into  Equation  (2.5)  yields; 


dipUi)  djpU.Uj)  dP  d 


dUi  dUj 


dx^  Bxj 


dxj  dXf 


Equation  (2.8)  is  also  known  as  the  “Boussinesq  Effective  Viscosity  Model”  since 
( |j,  +  ji, )  together  constitute  the  effective  laminar  and  turbulent  viscosities,  i.e., 


=  + 

As  mentioned  earlier,  the  turbulent  viscosity  varies  throughout  the  flowfield  depending 
on  the  local  flow  conditions.  It  represents  a  measure  of  both  eddy  size  and  strength,  and  must 
therefore  be  modeled  in  terms  of  parameters  which  reflect  these  two  quantities.  The  earliest 
attempt  to  develop  an  expression  for  turbulent  viscosity  that  varies  with  position  in  the 
turbulent  flow,  was  the  so-called  mixing-length  hypothesis.  This  was  first  initiated  by  Taylor 
(1915),  but  is  more  often  attributed  to  Prandtl  (1925),  some  10  years  later.  Prandtl,  who 
discovered  the  boundary  layer  in  1904,  now  reasoned,  partly  on  dimensional  grounds,  that  the 
turbulent  viscosity  was  proportional  to  the  magnitude  of  the  mean-strain  rate,  and  to  the  square 
of  a  characteristic  eddy-length  scale.  Since  the  mean  strain  rate  has  dimension  time'' ,  it  can  be 
thought  of  as  a  time  scale.  Therefore,  if  the  mixing  length  is  multiplied  by  this  quantity,  a 
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turbulent-velocity  scale  would  result.  The  expression  for  turbulent  viscosity  derived  from  the 
mixing-length  hypothesis,  and  written  below  as  Equation  (2.9),  represents  the  product  of  a 
length  scale  and  a  velocity  gradient.  The  eddy-viscosity,  according  to  the  mixing  length 
hypothesis,  is  then  prescribed  as  follows; 


[i,  =  p/^ 


dUi 

dxj 


(2.9) 


where,  /,  a  linear  dimension  called  the  mixing  length,  is  a  characteristic  of  the  flowfield,  and  is 
prescribed  empirically.  The  shortcomings  of  this  model  include  the  difficulty  in  prescribing  the 
mixing  length  for  complicated  flow  geometries  and  the  fact  that  jx,  vanishes  in  regions  of  zero 
mean  velocity  gradient.  Another  major  drawback  of  this  model  is  that  it  fails  to  take  into 
account  the  processes  associated  with  the  transport  of  turbulence,  thus  isolating  the  eddy- 
viscosity  at  one  particular  point  and  making  it  unresponsive  to  changes  in  the  flowfield  at  other 
locations.  In  other  words,  the  mixing  length  must  be  known  a  priori  everywhere  in  the 
flowfield  in  order  to  be  used.  Although  in  its  general  form,  the  mixing-length  hypothesis  is 
capable  of  producing  adequate  results  for  simple  boundaiy-layer-type  flows,  it  is  unable  to 
accurately  predict  flows  with  complex  strain  fields,  such  as  those  influenced  by  strong 
streamline  curvature.  Clearly,  for  complex  flowfields  such  as  those  encountered  behind  2-d  and 
3-d  geometries,  modifications  are  required  to  the  mixing-length  derivation  before  it  can  be  used 
successfully. 

Another  approach  to  the  turbulent  viscosity  dilemma  was  suggested  by  Kolmogorov  in 
1943  and  later  by  Prandtl  in  1945.  The  reasoning  is  that  there  exists  a  connection  between 
turbulence  parameters  and  the  eddy-viscosity.  They  formulated  a  revised  expression  for  the 
turbulent  viscosity,  known  as  the  Prandtl-Kolmogorov  turbulent-viscosity  formula  (1945); 
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(2.10) 


The  length  scale,  I,  is  now  proportional  to  the  size  of  the  turbulent  eddies,  and  must  be 
prescribed  explicitly.  The  turbulent  viscosity  is  considered  to  be  proportional  to  the  square  root 
of  the  TKE,  k,  which  may  be  calculated  throughout  the  flowfield,  and  easily  interpreted 
physically.  This  is  obtained  from  the  solution  of  a  related  transport  partial  differential  equation. 
A  conservation  equation  for  the  TKE  is  derived  by  multiplying  the  i"‘  momentum  equation  by 
M;  and  then  summing  from  i  =  1  to  3.  The  final  form  of  the  TKE  conservation  equation,  in  the 
fully  turbulent  regime,  for  flows  with  constant  density  in  tensor  notation  is  written  as  follows 
(see  also  Appendix  F); 


dk  -  dk  - dUi 

a7^^'ar=-“'“-ar 


ufu, 


pUj  dk  3m,  3m, 

+— -V  3— )-v(t^)(t^) 
p  dx.  dXj  dXj 


(2.11) 


1 


where  the  TKE  is  defined  as;  k  =  —  (MjM,  )  =  —  (m^  +  ),  and  v  =  |X  /  p 


The  general  terms  represented  in  Equation  (2.11)  describe  the  physical  processes  of 
turbulence.  The  first  term  on  the  right  is  the  production  term  and  represents  the  rate  at  which 
kinetic  energy  is  generated  and  transferred  to  the  turbulence.  The  second  term  on  the  right  is 
known  as  the  diffusion  term  and  contains  the  triple  velocity  correlation  terms  known  as  the 
turbulent  diffusion  transport  terms,  the  pressure  diffusion  terms  resulting  from  pressure  and 
velocity  correlations,  and  the  molecular  diffusion  terms.  The  third  term  on  the  right  is  the 
dissipation  term  and  describes  the  rate  at  which  turbulence  is  dissipated  into  heat  or  thermal 
energy.  The  first  term  on  the  left  describes  the  unsteady  behavior  of  TKE.  The  second  term  on 
the  left  is  the  convection  term  and  describes  the  rate  of  convection  of  TKE.  The  terms  on  the 
left  give  the  Eulerian  derivative  of  TKE  and  describe  its  rate  of  change  following  a  fluid 
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panicle.  For  closure  !o  lake  place,  Ihe  production,  dissipation,  turbulent  diffusion  (triple 
products),  and  pressure-velocity  correlation  terms  in  the  TKE  equation  have  to  be  modeled, 
since  they  involve  unknown  correlations.  This  closure  technique  requires  the  use  of 
approximations  which  are  for  the  most  part  deductions  and  assumptions  obtained  when 
studying  isolated  or  easily  defined  turbulent  flows. 

This  alternate  formulation  for  the  turbulent  viscosity,  together  with  similar  suggestions 
of  other  workers,  is  known  collectively  as  "one-equation  models."  They  are  now  distinguished 
from  the  models  discussed  previously,  which  are  often  referred  to  as  "zero-equation  models." 
Although  this  concept  overcomes  the  problem  of  the  eddy-viscosity  going  to  zero  as  the  mean 
velocity  gradient  vanishes,  an  a  priori  prescription  for  the  length  scale  of  turbulence,  /  ,  is  still 
necessary,  and  therefore  it  has  some  of  the  same  deficiencies  as  the  mixing-length  theory. 

To  reach  a  more  universal  model  that  overcomes  this  deficiency,  "two-equation 
models"  have  been  developed  where  the  transport  of  the  turbulence  length  scales  are  also 
accounted  for.  In  addition  to  the  Boussinesq  approximation,  and  the  TKE  equation,  a 
conservation  equation  for  a  variable  that  is  a  composite  of  k  and  I  is  solved.  This  constitutes  a 
complete  model  of  turbulence  where  no  prior  knowledge  of  turbulent  structure  is  necessary. 
The  differential  equation  for  an  independent  solution  for  I  is  very  awkward  and  not 
recommended.  Many  models  have  been  developed  in  the  past  that  used  various  possible 
combinations  of  k  and  / ,  such  as  the  /:  -co  model  of  Kolmogorov;  the  k-e  model  of  Chou 
(1945),  Davidou  (1961),  and  Harlow  and  Nakayama  (1968);  and  the  k-kl  model  of  Rotta 
(1968).  However,  the  most  popular  of  these  two-equation  models  is  the  k  —  E  model.  The 
standard  k  —  E  model  was  well  described  by  Jones  and  Launder  (1972),  but  widely  used  in  the 
form  of  Launder  and  Sharma  (1974)  who  “retuned”  the  closure  coefficients.  The  goal  is  to  find 
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and  solve  a  conservation  equation  for  the  dissipation  rate,  e,  with  suitable  closure 


approximations.  For  high  Reynolds-number  flows,  e  is  shown  to  be  related  to  the  turbulent 
length  scale  through; 


8  a 


/ 


(2.12) 


The  turbulence  closure  problem  is  now  resolved  and  the  need  to  specify  a  length  scale 
is  overcome.  The  turbulent  viscosity  can  then  be  obtained  by  dimensional  arguments  via; 


9. 


(2.13) 


However,  the  TKE  conservation  equation  and  the  much  more  complicated  dissipation 
rate  conservation  equation  contain  several  terms  that  are  practically  impossible  to  obtain,  such 
as  the  pressure-velocity  correlations  terms,  and  the  double  and  triple  correlations  of  fluctuating 
velocity  terms.  These  terms  are  generally  labeled  in  a  fashion  similar  to  the  TKE  terms 
mentioned  earlier,  but  in  this  case  instead  of  k ,  they  identify  e.  Therefore,  in  order  for  closure 
to  be  possible,  certain  simplifying  assumptions  about  these  terms  have  to  be  made.  For 
instance,  measurements  made  in  simple  turbulent  flows  give  the  three  “universal”  constants. 
Cel,  Ce2,  and  Cji .  Ce2  is  obtained  from  measurements  of  decaying  turbulence  similar  to  that 
found  behind  a  grid  in  the  absence  of  mean  velocity  gradients.  Also,  from  experimental 
observations  made  in  the  near-wall  region  of  a  fully  developed  turbulent  boundary  layer,  both 
convection  and  diffusion  of  TKE  are  very  small  and  may  be  neglected,  which  means  that 
turbulence  generation  is  relatively  in  balance  with  turbulence  dissipation.  Under  these 
conditions,  the  TKE  for  eddy  viscosity  reduces  to  the  PrandtTs  mixing-length  formulation. 
This  allows  the  approximation  of  the  other  two  constants,  Cei,  and  . 
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The  k-e  turbulence  model  has  proven  to  be  useful  for  a  large  variety  of  flows 
including  those  where  the  streamlines  do  not  curve  significantly  in  shear  regions,  such  as 
straight  boundary  layers,  wakes  and  jets,  where  the  shear-stresses  dominate  the  coupling 
between  the  turbulence  and  mean-velocity  fields.  However,  in  curved  flows,  particularly  when 
the  normal-stresses  contribute  greatly  to  mean  transport,  such  as  in  recirculating  flow,  the 
k-e  model  fares  less  well.  Although  the  general  features  of  the  flow  can  sometimes  be 
captured  reasonably,  the  predicted  turbulence  quantities  often  show  poor  agreement  with 
experimental  measurements. 

One  of  the  major  drawbacks  of  the  two-equation  turbulence  models  and  for  that  matter 
any  eddy-viscosity-based  turbulence  model  is  the  implicit  assumption  of  isotropic  turbulence. 
That  is,  the  models  do  not  account  for  rotational  or  translational  differences  in  the  turbulent 
normal  stresses.  This  is  due  to  the  turbulent  pressure  term  which  tends  to  overwhelm  the 
velocity  gradient  terms  in  most  flows  and  hence  forces  the  three  normal  stresses  to  2k/3.  In 
algebraic  models,  the  turbulent  viscosity  is  a  local  function  while  in  the  two-equation  models  it 
is  a  more  complicated  and  more  general  function  governed  by  two  partial  differential  equations. 
The  eddy-viscosity  model  simply  relates  the  Reynolds  stresses  and  the  local  mean  velocity 
gradients  or  rates  of  mean  strain  via  some  empirically  determined  turbulent  viscosity 
relationship. 

It  is  well  known  that  all  shear  flows  and  wall  turbulence  flows  are  anisotropic  which 
raises  doubts  in  the  use  of  eddy-viscosity  based  turbulence  models.  To  remedy  this  problem, 
Reynolds  Stress  Models  (RSM)  have  been  proposed.  These  models  complete  closure  to  the 
Reynolds  equations  without  using  the  turbulent  viscosity  assumption.  Like  the  two-equation 
models,  the  shortcoming  of  this  modeling  technique  is  the  need  to  make  assumptions 
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concerning  the  various  terms  to  be  evaluated,  such  as  those  encountered  in  the  TKE  transport 
equation.  The  terms  that  have  been  measured  in  this  study  include  the  nine  third-order  turbulent 
velocity  correlations. 

All  the  models  mentioned  above  rely  on  assumptions  and  approximations  for 
evaluation  of  high  order  correlations.  Their  success  remains  closely  related  to  how  well  these 
terms  are  modeled  and  whether  they  are  important  in  determining  the  transport  of  turbulence. 
Clearly,  there  still  exists  plenty  of  room  for  improvement  and  refinement. 

Among  the  non-eddy  viscosity  based  numerical  computational  methods  available  is 
the  Algebraic  Stress  Model  (ASM).  This  technique,  first  proposed  by  Rodi  (1972),  is 
sometimes  considered  simple,  and  yet  still  efficient  to  predict  complex  turbulent  shear  flows. 
With  this  method,  the  differential  transport  equations  for  the  Reynolds  stresses  are  replaced 
by  a  set  of  algebraic  equations  for  Reynolds  stresses  assuming  that  the  convection-diffusion 
term  of  the  kinetic  energy  equation  is  proportional  to  that  of  the  Reynolds-stress  transport 
equations. 

Of  particular  interest  is  the  Renormalization  Group  (RNG)  theory  first  applied  to  the 
k-e  turbulence  model  by  Yakhot  and  Orszag  in  1986,  and  later  improved  by  Yakhot,  et  al.  in 
1992.  The  goal  of  the  RNG  theory  is  to  systematically  remove  the  smallest  turbulent  scales 
to  a  point  where  the  remaining  scales  are  resolvable.  The  constants  used  in  the  RNG  k-£ 
model  are  derived  from  theory,  removing  much  of  the  ad-hoc  nature  of  the  standard  k-e 
model,  and  represents  the  major  functional  difference  between  the  models.  Also,  unlike  the 
Reynolds  Stress  models  which  are  based  on  Reynolds  averaging,  a  technique  that  smoothes 
out  some  of  the  important  features  of  turbulence,  the  RNG  theory  systematically  removes  the 
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.smulle,sl  scales  of  lurbuicnce  lo  (he  pom(  where  the  remaining  scales  am  .esolvable  with 
currently  available  computer  capabilities  (Yakhot  and  Orszag,  1992). 

Of  all  modeling  techniques  mentioned  above,  the  twoequation  eddy  viscosity 

turbulence  model  remains  the  most  widely  used  due  to  its  simplicity  and  flexibility.  Although  it 

does  not  properly  model  the  turbulence  parameteis  in  highly  complex  flows,  it  gives  results 

good  enough  to  allow  the  engineer  to  make  certain  decisions.  Some  recently  modified  veisions 

of  the  k  -  e  model  have  given  better  results  than  the  standard  model  as  will  be  shown  later  in 
this  study. 

The  current  study  shows  a  comparison  between  the  experimental  data  and  results 
from  three  numerical  solutions;  1)  a  symmetric,  steady-state  k-e  solution;  2)  a  standard  k-e, 
time-accurate  solution;  and  3)  an  RNG  k-e,  time-accurate  solution.  All  solutions  are  two 

dimensional,  and  the  mean  velocities  and  Reynolds  stresses  are  compared  to  the 
experimental  measurements. 
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CHAPTER  3 


FYPERTMENTAT.  APPARATUS 


3.1  Facility 

Schematics  of  the  wind  tunnel  and  the  2-d  test  section  used  for  this  study  are  shown 
in  Figure  1.  The  test  facility  is  located  in  Test  Cell  18  at  the  Aero  Propulsion  and  Power 
Directorate  of  Wright  Laboratory,  Ohio,  and  was  designed  and  constructed  to  accommodate 
flow  visualization  and  laser-based  diagnostics.  Full-view  fused  quartz  windows  on  the  top, 
bottom,  and  the  two  side  walls,  provide  optical  access  to  the  76.2-mm  x  127-mm  x  762-mm 
(3-in  X  5-in  x  30-in)  test  section.  Air  flow  was  supplied  by  a  7.5-horsepower  radial  fan 
blower  followed  by  a  152.4-mm  (6-in)  diameter  honeycomb  flow  straightening  section  and 
settling  chamber.  The  flow  straightener  was  placed  close  to  the  fan  outlet  to  reduce  the  large 
scale  turbulence  produced  by  the  fan.  The  air  was  further  conditioned  by  passing  through 
another  set  of  honeycomb  flow  straighteners  placed  304.8  mm  (12  in)  upstream  of  the  test 
section  entrance.  Fine  mesh  screens  were  removed  from  this  flow  downstream  of  the 
straightening  section  because  these  became  clogged  with  seed  particles  in  an  unpredictable 
and  nonuniform  fashion.  Test  section  flow  parameters,  such  as  static  pressure,  and  total 
temperature  were  monitored  and  recorded  by  a  micro-computer  based  data  acquisition 
system.  For  the  present  study,  the  inlet  centerline  reference  velocity,  was  monitored 

using  a  pitot/static  tube,  connected  to  a  high  resolution  manometer,  and  located  at  xfll=-24 
(12  inches  upstream  of  the  bluff  body).  The  averaged  inlet  reference  velocity  was  -23.78  m/s 
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(78  ft/s),  corresponding  to  a  tunnel  Mach  number  of  0.08  and  a  Reynolds  number  of  2.0  x 
10^  based  on  inlet  flow  conditions  and  the  bluff-body  height,  H. 


b)  Test  Section 


Figure  1.  Experimental  test  facility. 

3.2  Bluff  Body 

The  shape  of  the  bluff-body  flameholder  affects  the  flow  stability  characteristics 
through  the  influence  on  the  size  and  shape  of  the  wake  region  (Lefebvre,  1983).  One  of  the 
pnmary  functions  of  a  flameholder  is  its  ability  to  anchor  the  flame  by  creating  a 
recirculation  zone  that  enhances  the  mixing  of  the  injected  fuel  and  the  incoming  air.  The 
recirculation  zone  is  created  when  the  static  pressure  in  the  central  core,  just  downstream  of 
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the  bluff  body,  becomes  low  enough  to  create  flow  recirculation.  The  siz»  of  flie  recirculation 
zone  can  be  increased  by  modifying  the  flameholder  geometry.  Increasing  the  lecirculation 
zone  length  has  a  tendency  to  raise  the  maximum  reverse  mass  flow  and  tiie  rate  of 
entrainment,  albeit  with  some  increase  in  pressure  loss  (Lefebvre,  1983).  According  to  Fox 
&  McDonald  (1978),  a  hemispherical  geometry,  with  the  rounded  head  facing  upstream, 
showed  the  lowest  total  (pressure  and  friction)  drag  coefficient,  typically  a  of  0.38. 
Therefore,  a  bluff  body  with  a  blunt  nose  geometry  was  selected  for  this  study.  A  schematic 
of  the  bluff-body  model  and  its  orientation  to  the  coordinate  system  is  shown  in  Figure  2. 
The  25.4-mm  (1-in.)  long  model  consisted  of:  1)  a  6.35-mm  (0.25-in)  radius  hemispherical 
nose,  2)  a  15.875-mm  (0.625-in)  thick  rectangular  center  section  and,  3)  a  3.175-mm  (0.125- 
in)  base/dump  plate.  This  last  section  is  interchangeable  to  allow  for  attachment  of  various 
vortex  generator  designs.  The  bluff-body  height  (H)  was  12.7  nun  (0.5  in),  and  spanned  the 
entire  test  section  giving  a  tunnel  blockage  of  16%. 


~T 

0.5  " 

A. 

a)  Isometric  View  b)  Side  View 

Figure  2.  Bluff-body  model  schematic. 

The  model  was  mounted  directly  to  the  side  wall  windows  as  indicated  in  Figure  lb, 
in  order  to  eliminate  flow  disturbing  model  supports  (i.e.  stings).  The  model  was  mounted 
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127  mm  (4  inches)  into  the  test  section,  downstream  of  the  side  window  leading  edge.  Great 
care  was  taken  to  minimize  the  effect  of  the  upstream  influences  on  downstream  flow 
behavior,  since  previous  work  had  shown  that  the  flow  in  the  wake  of  a  bluff  body  can  be 
three-dimensional.  This  arrangement  (see  Figure  lb)  provided  the  best  optical  access  for 
probing  the  entire  flowfield  surrounding  the  bluff  body.  Care  was  taken  to  align  the  bluff 
body  with  respect  to  the  laser  Doppler  velocimeter  probe  volume  by  running  the  laser  beam 
along  the  edge  of  the  bluff  body. 


3.3  Laser  Doppler  Velocimetrv 

Laser  Doppler  velocimetry  (LDV)  has  been  used  for  fluid  flow  measurement  since 
the  mid  1960’s.  Durst,  et  al.  (1976)  describes  the  theory  and  application  of  LDV  at  length. 
LDV  is  a  nonintrusive  instrument  used  to  measure  turbulent  velocity  fields  without 
disturbing  the  flow,  and  where  conventional  techniques  (i.e.,  hot  wire)  are  not  practical. 
Typical  applications  include  turbulent  boundary  layers,  vortices  and  wakes,  such  as  the 
recirculation  zone  behind  a  bluff  body  or  a  backward-facing  step.  A  single  probe  volume  is 
created  by  crossing  two  polarized  laser  beams.  At  the  crossing,  the  two  beams  constructively 
and  destructively  interfere  to  create  fringes  of  light  as  shown  in  Figure  3.  As  seed  particles  in 
the  flow  pass  through  the  probe  volume,  the  scattered  light  from  the  fringes  yield  an 
amplitude  modulated  sinusoidal  signal,  called  a  Doppler  burst.  The  frequency  of  the 
scattered  light  as  the  particle  crosses  bright  and  dark  fringes  is  related  to  the  particle  velocity 
component  perpendicular  to  the  fringe  direction  by  the  following  relation; 


U{x,y,z)  = 


2sin((|)  /2) 


(3.1) 
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Where  X  is  the  laser  beam  wavelength,  the  measured  instantaneous  Doppler  frequency  of 
the  particle,  and  <|)  is  the  angle  between  the  two  laser  beams  being  crossed  to  form  the  probe 
volume,  also  identified  as  the  beam  separation  angle  (see  Appendix  A). 


Figures.  Typical  dual  beam  laser  Doppler  anemometer 

It  is  important  to  note  that  the  velocity  measured  is  the  component  perpendicular  to 
the  fringes  in  the  probe  volume.  Two-  and  three-dimensional  LDV  systems  are  created  by 
focusing  additional  pairs  of  laser  beams  of  different  colors  at  the  same  location,  but  at 
another  orientation.  By  varying  the  orientation  of  the  three  sets  of  beams,  three  independent 
components  of  the  same  velocity  vector  can  be  measured  simultaneously.  The  scattered  light 
is  collected  by  a  lens  and  focused  onto  a  pinhole  aperture  of  a  photomultiplier  tube  (PMT) 
(or  photodetector),  which  transmits  an  electrical  signal;  a  voltage  modulated  at  the  Doppler 
frequency  that  can  be  analyzed  by  a  signal  processor.  The  LDV  signal  processors  amplify 
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and  filter  the  signals  from  the  photomultiplier,  validate  the  Doppler  frequency  samples,  and 
finally  compute  the  Doppler  period  which  is  the  reciprocal  of  the  Doppler  frequency.  This 
integer  number  (in  nanoseconds)  is  then  transmitted  to  the  computer.  The  Doppler  frequency 
of  the  signal  gives  the  velocity  of  the  particle  at  the  instance  it  passes  through  the  probe 
volume.  The  particles  arrive  randomly  at  the  measurement  volume,  and  hence  the  output 
signal  of  an  LDV  used  in  a  turbulent  flow  is  a  randomly-sampled  signal. 

Both  one-component  and  three-component  LDV  system  configurations  were  used  in 
this  study.  The  one-component  system  was  used  for  validation  of  the  three-component 
system  as  well  as  for  mass-balance  verification.  The  orientation  of  this  system  was  such  that 
the  beams  passed  through  the  side  windows  of  the  test  section,  parallel  to  the  top  and  bottom 
tunnel  surfaces,  and  thus  allowed  for  the  measurement  of  the  axial  velocities  at  locations 
very  close  to  the  top  and  bottom  tunnel  walls.  The  far-field  measurements  were  taken  using 
two  separate  TSI  fiber-optic  probes  based  LDV  systems.  The  three-component  LDV  system 
was  mounted  around  the  rectangular  test  section  and  was  operated  in  the  forward-scatter 
mode.  Simultaneous  measurements  of  axial  (U),  transverse  (V),  and  spanwise  (W) 
components  of  the  velocity  vector  (see  Appendix  B)  were  made  throughout  the  tunnel  test 
section.  In  both  the  near-  and  the  far-field  measurements,  the  third  component  came  from  a 
second  independent  one  Watt  argon-ion  laser.  The  optics  were  arranged  to  obtain  a  minimum 
beam  waist  diameter  at  the  probe  volume  location.  This  gave  the  highest  beam  intensity.  The 
transmitting  and  receiving  optics  were  mounted  on  an  xyz  optical  table  capable  of  placing 
the  measurement  probe  volume  at  various  points  throughout  the  tunnel  test  section  with  a 


resolution  of  ±0.03  mm. 


3.3.1  Far-Field  LDV  System 


A  three-component  LDV  system  operating  in  the  forward-scatter  mode  was  mounted 
around  the  rectangular  test  section  to  make  simultaneous  measurements  of  the  U,  V,  and  W- 
velocity  components.  The  optical  system  used  in  the  experiments  involved  two  separate 
fiber-optic  based  TSI  Inc.  laser  Doppler  systems.  This  allowed  the  laser  beams  to  be  oriented 
orthogonal  to  each  other  giving  direct  measurements  of  the  three-orthogonal  velocity 
components.  A  TSI  Colorburst  beam  separator  was  used  to  separate  the  various  wavelengths. 
The  two  most  intense  argon  laser  lines,  the  green  (A.=514.5  nanometers)  and  the  blue  (A.=488 
nanometers),  were  split  into  two  beams  of  near  equal  intensity,  frequency  shifted,  and 
directed  via  fiber-optic  cables  to  a  TSI  model  9277  two-component  fiber-optic  probe.  The 
transmitting  optics  included  a  3.75x  beam  expander  (for  a  full  explanation  of  this  term  and 
other  laser  related  terms,  see  Appendix  C)  and  had  a  focal  length  of  450  mm.  The  orientation 
of  this  system  was  such  that  the  beams  passed  through  the  side  windows  of  the  test  section 
and  thus  allowed  for  the  measurement  of  the  axial  and  transverse  (U,V)  velocity 
components.  In  the  other  LDV  system,  the  violet  argon  laser  line  (X=476.5  nanometers)  from 
a  second  argon  laser  was  directed  via  fiber-optic  cable  to  a  TSI  model  9274  one-component 
fiber-optic  probe.  This  allowed  the  measurement  of  the  spanwise  velocity  component  (W). 
The  orientation  of  this  system  was  such  that  the  beams  were  orthogonal  to  each  other  and 
thus  passed  through  the  top  and  bottom  test  section  windows.  In  order  to  avoid  directional 
ambiguity  that  would  otherwise  result  from  the  highly  turbulent  recirculating  flow,  a  40- 
MHz  frequency  shift  was  applied  to  one  beam  in  each  of  the  three  components.  Both  argon- 
ion  lasers  were  operated  at  approximately  one  Watt  for  all  measurements.  The  approximate 
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probe  volume  dimensions  were  600-|xm  long  and  SO-pm  diameter  for  the  two-component 
system  and  l5()-pm  long  and  l(X)-pm  diameter  for  the  one-component  system  (see  Appendix 
A).  The  fringe  spacings  for  the  green,  the  blue,  and  the  violet  were  1.738-|Xm,  1.827-pm,  and 
3.439-pm,  respectively.  The  receiving  optics  consisted  of  a  120-nim  lens  focused  onto  a 
0.25-mm  diameter  photomultiplier-tube  aperture  via  an  appropriate  dichroic  filter  to 
selectively  pass  either  the  blue  or  green  beams.  Additionally,  laser  line  filters  were  mounted 
to  each  photomultiplier  tube.  The  violet  signal  was  collected  in  the  forward-scattered  mode 
which  was  at  ninety  degrees  with  respect  to  the  blue  and  the  green  axis  in  this  system.  A 
476.5-nm  laser  line  filter  was  used  here  to  minimize  signal  cross-talk.  Due  to  the 
requirement  imposed  by  the  processing  electronics,  that  axial,  transverse,  and  spanwise 
velocity  components  be  obtained  simultaneously,  a  great  effort  was  made  to  ensure  that  the 
three  separate  probe  volumes  were  coincident  by  passing  all  three  through  a  35-pm  pinhole 
placed  at  45  degrees  to  the  violet  and  blue/green  beam  propagation  directions.  The  effective 
measurement  probe  volume  shape  now  becomes  a  sphere  with  a  diameter  equal  to  the  beam 
waist  with  a  cross-sectional  area  approximately  equal  to  0.008  nmi^,  that  of  the  violet  beam. 
Although  the  effective  probe  volume  size  is  reduced  by  the  addition  of  the  third  component, 
the  set-up  had  to  be  perfected  in  order  to  ensure  that:  1.)  all  six  beams  passed  through  the 
pinhole,  and  2.)  that  the  beam’s  intensity  was  optimized  as  well  as  other  factors  which  have 
to  do  with  optimizing  the  data  rate.  These  efforts  are  all  responsible  for  producing  the  high 
quality,  smooth  and  symmetric  profiles  that  were  observed  in  the  resulting  measurements  of 
this  study. 


3.3.2  Near-Field  LDV  System 


The  system  described  above  has  a  severe  limitation  insofar  as  it  is  geometrically 
impossible  to  obtain  data  in  the  near-wake  region.  It  is  physically  impossible  to  get  all  three 
probe  volume  to  coincide  at  a  distance  less  than  12.7  mm  (0.50  in)  from  the  base  of  the  bluff 
body,  when  at  the  center  span  location.  Therefore,  a  different  system  .setup  was  used  in  the 
near-field  region.  In  this  setup,  the  green  and  blue  fringes  were  oriented  at  90  degrees  with 
respect  to  each  other  (i.e.  orthogonal)  and  oriented  at  45  and  135  degrees  with  respect  to  the 
horizontal  axis  of  the  test  section.  This  allowed  the  measurements  to  be  made  to  within  6.35 
mm  (0.25  in)  of  the  trailing  edge  of  the  bluff  body.  The  transformation  equations  discussed 
in  Appendix  B  were  used  to  recover  the  velocity  components.  The  LDV  used  for  this  set  of 
experiments  was  a  conventional  two-color,  four-beam,  back-scatter  system  (as  shown  in 
Figure  4).  However,  to  achieve  higher  data  rates  and  to  increase  the  signal-to-noise  ratio,  the 
receiving  optics  were  arranged  in  a  forward-scatter  configuration  for  this  study.  The  LDV 
system  included  Bragg-cell  modulators  to  allow  a  net  frequency  shift  of  40  MHz  in  both  the 
green  and  the  blue  beams.  No  beam  expander  was  used  in  this  system  set-up.  An  80-mm 
diameter,  181.5-mm  focal  length  lens  for  transmitting,  was  used.  The  blue  and  the  green 
lines  of  an  argon-ion  laser  operating  at  one  Watt  were  used  in  all  measurements.  The 
approximate  measurement  volume  dimensions  based  on  1/e^  intensity  points  of  the 

measuring  volume  were  600  |im  in  length  and  80  |im  in  diameter  (see  Appendix  A). 
Scattered  light  from  the  probe  volume  was  separated  with  a  dichroic  filter  so  that 
approximately  80%  of  each  color  was  passed  to  its  respective  photomultiplier  tube.  Narrow 
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biindpass  filters  arc  used  to  Further  reduce  signal  cross  talk.  The  violet  laser  line  was  brought- 
in  from  the  top  similar  to  the  far-field  measurements  case  discussed  earlier. 


For  both  LDV  systems,  in  order  to  ensure  spatial  coincidence  of  the  three  probe 
volumes,  a  35-^m  pinhole  was  placed  at  45-degrees  at  the  focal  point  and  all  six  beams  were 
accurately  co-located  by  forcing  (steering)  them  to  pass  through  the  pinhole.  The  entire 
optical  system  was  mounted  on  a  three-axis  traversing  table,  which  allowed  the  probe 
volume  to  be  positioned  with  a  resolution  of  0.03  mm  (0.001  in).  The  origin  (x,y,z=0)  for  all 
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LDV  measurements  was  located  on  the  transverse  and  spmwise  centerline  at  the  bailing 
edge  of  the  flameholder  as  shown  in  Figures  1  and  2.  Each  LDV  channel  measured  a  yelocity 
component  for  particles  traveling  through  the  probe  volume. 

3.4  Flow  Seeders 

According  to  Jones,  et  al.  (1990),  seeding  is  the  most  difficult  and  often  the  most 
neglected  aspect  of  laser  velocimetry.  Since  particles  in  the  flow  are  the  source  of  the 
Doppler  signal,  seeding  becomes  an  important  factor  in  the  performance  of  an  LDV  system. 
A  classical  dilemma  arises  due  to  the  necessity  to  provide  particles  large  enough  to  scatter 
sufficient  light  for  good  signal  quality  so  as  to  be  detected  by  the  photodetector,  yet  small 
and  light  enough  to  closely  follow  curved  flow  streamlines  even  in  fluid  flow  regimes  having 
steep  velocity  gradients.  Small  seeding  particles  are  also  necessary  in  order  to  capture 
faithfully  the  turbulent  fluctuations,  and  the  higher  order  velocity  correlations  (Maurice, 
1992).  Any  seeding  solution  represents  a  trade-off  that  must  be  made  between  signal  quality 
and  particle  motion  fidelity. 

Two  different  flow  seeders  were  used  here,  one  for  the  far-field  and  one  for  the  near¬ 
field  measurements.  In  both  cases,  the  seed  particles  were  injected  upstream  of  the  settling 
chamber  in  order  to  ensure  a  uniform  seed  density  at  the  inlet  to  the  test  section.  This  also 
provided  seed  densities  sufficient  to  give  data  validation  rates  in  excess  of  10,000  per  second 
for  each  channel  (green,  blue,  violet  beams)  over  most  of  the  flowfield.  The  coincidence  data 
rate,  for  both  studies,  varied  between  1000/sec  in  the  recirculation  zone  and  5000/sec  in  the 
freestream  and  the  far-wake  region  of  the  bluff  body. 
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3.4.1  Far-Field  Seeding 

For  the  far-field  measurements,  a  chemical  seeder  developed  by  Craig,  et  al.  (1984) 
which  produced  -1-jxm  diameter  titanium  dioxide  particles,  was  used.  The  titanium  dioxide 
(TiO^)  seed  particles  were  generated  by  reacting  titanium  tetrachloride  {TiCl^)  with  moist 
air,  according  to  the  following  chemical  reaction; 

TiCl,  +  ->  TiO^  +  A{HCl)  (3.2) 

The  titanium  dioxide  seed  particles  produced  here  are  the  same  as  those  commonly 
found  in  “toothpaste.”  However,  the  other  product  of  this  chemical  reaction,  is  Hydrochloric 
acid,  and  is  well  known  to  be  a  health-hazard.  Therefore,  care  was  taken  to  ensure  that  small 
quantities  of  seed  particles  were  generated  and  properly  disposed  of  by  venting  the  products 
through  a  roof-vent.  The  quantity  was  measured  and  controlled  via  an  in-line  rotometer. 
Titanium  tetrachloride  must  be  extremely  well  contained,  because  of  the  spontaneous 
reaction  that  takes  place  when  placed  in  contact  with  moisture  in  the  atmosphere. 

The  reaction  chamber  was  designed  using  the  specifications  given  by  Nejad  (1986). 
Craig,  et  al.  (1986)  measured  the  particle  sizes  generated  by  this  device  and  found  that  they 
were  fairly  spherical  and  in  the  0.2  -1  |xm  diameter  range.  The  reaction  chamber  was  capable 
of  supplying  seed  densities  sufficient  to  give  data  validation  rates  -5000  per  second.  Also, 
the  (Ti02)  particles  coated  the  test  section  windows  very  quickly,  especially  in  the 
recirculation  zone,  and  consequently  reduced  the  data  validation  rate.  The  validation  data 
rate  dropped  even  further  when  the  measurements  were  taken  in  the  recirculation  zone, 
where  fewer  seed  particles  are  found  and  turbulence  intensities  are  higher. 
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3.4.2  Near-Field  Seeding 


For  the  near-field  measurements,  a  particle  generator  designed  by  Rabe  and  Sabroske 
(1994)  was  used.  The  reference  gives  a  clear  description  of  the  seeder.  The  generator  uses 
four  Laskin  nozzles  submerged  in  liquid  glycerin.  Operational  guidelines  given  in  their  work 
were  followed.  They  measured  the  particle  size  distribution  and  found  that  they  ranged  from 
0.5  to  1.0  p.m  in  diameter. 

In  both  the  far-  and  near-field  measurements,  data  acquisition  was  stopped  and  the 
test-section  windows  were  cleaned  when  the  data  rate  dropped  below  1000  per  second.  Away 
from  the  bluff  body  and  in  the  freestream,  the  coincident  data  rate  appeared  to  remain ' 
relatively  unchanged,  at  5000  per  second. 

The  smooth  measured  profiles  indicate  that  the  sample  size  of  5120  was  sufficient  to 
obtain  good  statistical  estimates  for  the  mean  and  higher  order  moments  and  that  the  flow 
was  statistically  steady.  The  smoothness  was  further  demonstrated  when  the  data  were 
numerically  differentiated  in  the  axial  and  the  transverse  directions,  using  central  ‘ 
differencing  technique.  The  differentiated  data  was  used  in  the  turbulent  kinetic  energy 
equation,  and  the  transport  equation. 

3.5  Data  Acquisition  and  Proce.s.sing 

The  data  acquisition  and  processing  system  consisted  of  a  three-channel  IFA-750 
Digital  Burst  Correlator  and  Analyzer  by  TSI,  Inc.  The  IFA-750  signal  processor  allows  for 
measurement  of  the  Doppler  frequencies  of  all  three  velocity  components  and  ensures  that 
the  measurements  are  coincident  to  within  a  user-specified  coincidence  window.  For 
multichannel  operation,  the  signal  processor  ensures  that  all  active  channels  produce  a  valid 
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data  point  within  the  coincidence  window  before  data  are  accepted.  Once  the  point  is 
accepted,  data  are  stored  and  then  multiplexed  into  the  computer  via  a  data  parallel  interface. 

The  IFA-750  uses  the  signal-to-noise  ratio  (SNR)  to  distinguish  the  signal  burst  from 
the  noise.  It  receives  electrical  signals  from  the  photomultipliers,  and  extracts  Doppler 
frequency  information.  The  data  are  then  analyzed  using  the  Flow  Information  Display 
(FIND)  software  package  which  displays  the  turbulence  statistics  in  real-time. 

The  Doppler  signals  produced  by  an  LDV  occur  in  short  bursts,  at  random  time 
intervals  with  typically  10  to  50  cycles  in  each  burst.  They  are  processed  after  passing 
through  preset  low  and  high-pass  filters  on  each  channel  processor.  There  are  three  types  of 
signals:  1.)  A  true  Doppler  signal  whose  frequency  is  proportional  to  the  fluid  velocity  being 
measured,  2.)  a  pedestal,  and  3.)  wide-band  noise. 

The  persistent  wide-band  background  noise,  the  pedestal,  and  the  fact  that  the 
Doppler  signals  occur  at  random  intervals,  all  make  processing  the  real  signal  a  “challenging 
task”.  In  order  to  extract  the  Doppler  frequency,  the  Doppler  signal  has  to  be  separated  from 
the  pedestal  and  the  background  noise  has  to  be  reduced  as  much  as  possible.  This  is  done 
through  the  high-pass  filter  which  eliminates  the  pedestal,  while  the  low-pass  filter  reduces 
the  amplitude  of  the  wide-band  or  the  background  noise.  Also,  the  signals  occurring  at 
random  intervals  in  short  bursts  are  detected  by  using  a  Burst  Detector  and  Sampler.  The 
burst  detector  is  used  to  discriminate  between  the  Doppler  signal  and  the  noise.  When  a 
signal  exceeds  a  preset  (SNR)  level,  a  burst  timer  measures  the  length  of  the  signal  burst. 
Other  steps  follow  to  ensure  that  only  good  measurements  are  validated  and  sent  to  the 
computer.  A  more  detailed  discussion  of  how  the  IFA-750  extracts  velocity  information  from 
Doppler  signals,  is  available  through  publications  such  as  the  TSI  manual  (1992). 
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The  TSI  IFA-750  autocorrelator  was  used  to  process  the  LDV  data,  and  extract 
Doppler  frequencies.  The  system  was  set  to  operate  in  the  coincidence  mode,  validate  only  a 
single  measurement  per  burst,  and  to  count  at  least  eight  fringes.  The  coincidence  window 
was  set  to  20  psec  throughout  this  study.  This  means  that  coincidence  between  measured 
velocity  components  (U,  V,  and  W)  has  been  defined  to  take  place  when  all  three  laser 
velocimeter  channels  send  a  validated  velocity  measurement  to  the  signal  processor  buffer 
within  20  psec  of  one  another.  It  is  considerably  more  difficult  to  obtain  physically 
meaningful  three-component  data  than  two-component  data.  This  is  largely  due  to  the  difficulty 
in  achieving  and  maintaining  coincidence  of  the  three  probe  volumes  to  obtain  coincident 
three-component  data.  This  difficulty  caused  a  drop  in  the  coincidence  data  rate  by  a  factor  of  3 
to  9  less  than  the  lowest  of  the  three  individual  data  rates  was  noted.  The  low  and  high-pass 
filters  were  set  to  20  MHz  and  100  MHz,  respectively,  on  each  of  the  three  processors. 

In  the  post-processing  of  an  LDV  data  set,  “bad”  data  points  may  be  encountered  and 
should  be  rejected.  No  matter  how  sophisticated  the  electronic  circuitry  of  the  autocorrelator 
is,  there  may  be  some  validated  points  which  are  not  “true”  particle  velocities,  but  rather  a 
result  of  signal  noise  such  as  optical,  electronic,  or  a  combination  of  both.  These  points  are 
called  “outliers”  or  “spurious”  since  they  often  lie  many  standard  deviations  from  the  mean. 
Such  points,  though  few  in  number,  can  have  a  significant  effect  on  the  magnitude  of  the 
calculated  standard  of  deviation  of  the  data  set.  Therefore,  in  this  study,  extreme  outliers 
were  removed  first.  Then  the  mean  and  standard  deviation,  a,  were  calculated.  All  data  lying 
outside  a  designated  multiple  of  o  were  then  eliminated  from  the  original  data  set.  After  that, 
the  remaining  data  are  considered  as  the  ‘true’  data  set,  for  which  a  new  mean  and  a  new 
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standard  deviation  are  calculated.  In  addition  to  the  frequency  data,  the  time  between  data 
and  the  probe  volume  particle  residence  time  were  recorded.  All  sample  sets  included  5120 
individual  realizations  for  each  velocity  component.  The  data  were  filtered  using  a  ±3a  cut¬ 
off  limit  before  final  statistical  parameters  were  calculated.  In  most  cases,  the  signal-to-noise 
ratio  was  high  and  few  samples  were  discarded.  The  statistical  uncertainties  of  all  measured 
parameters  were  also  calculated  using  the  jackknife  method  (see  Benedict  &  Gould,  1996). 

3.6  Velocity  Bias 

The  problem  of  velocity  bias  in  LDV  measurements  has  been  thoroughly  examined 
since  it  was  first  recognized  by  McLaughlin  and  Tiederman  (1973).  Velocity  bias  arises  from 
the  fact  that  signal  processors  make  discrete  velocity  measurements  from  individual 
realizations  of  seed  particles  passing  through  the  probe  volume.  In  a  uniformly  seeded  flow, 
the  number  of  particles  per  unit  time  passing  through  the  measurement  volume  is 
proportional  to  the  flow  rate  through  that  volume.  Therefore,  simple  arithmetic  averaging  of 
an  ensemble  of  realizations  will  produce  a  mean  velocity  biased  towards  values  greater  than 
the  true  temporal  mean  velocity.  This  effect  is  most  pronounced  in  highly  turbulent  flows, 
and  a  number  of  different  schemes  have  been  proposed  for  its  removal.  No  velocity  bias 
corrections  were  used  in  the  present  study. 

3.7  Data  Reduction  Packa£e 

An  Intel  486-based  PC  data  analysis  system  was  used  for  data  acquisition,  evaluation, 
manipulation,  and  storage.  In  order  to  collect  a  realistic  statistical  sample,  a  minimum  of 
5120  realizations  per  channel  were  taken  at  each  measurement  location.  The  autocorrelator 
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time-stamps  each  velocity  realization  to  allow  for  post-processing  of  the  data.  The  data 
analysis  software  supplied  by  TSI,  Flow  Information  Display  (FIND),  was  found  to  be 
inadequate  for  this  study  for  the  following  reasons;  1.)  it  did  not  allow  control  over  angles 
between  laser  beams  for  the  various  components,  2.)  it  did  not  allow  control  over  the  Bragg 
shift  direction,  3.)  a  point-to-point  normalization  of  data  was  needed  to  accommodate  the 
variations  in  the  reference  velocity,  and  4.)  it  did  not  calculate  statistical  uncertainties. 

In  order  to  overcome  these  problems,  a  new  data  reduction  computer  program 
developed  by  Gould  (1995),  was  used.  This  specialized  software  reads  the  standard  TSI-750 
generated  output  file  that  contains  header  information  and  the  velocity  and  time  information  for 
each  recorded  velocity.  User  supplied  information  for  the  total  temperature,  pressure, 
differential  pressure,  and  reference  velocity  at  each  point  (x,  y,  z)  are  then  read  from  an  outside 
file.  Advanced  LDV  data  reduction  techniques  are  then  used  to  process,  analyze,  and  prepare 
the  data  into  a  presentable  format  (a  detailed  derivation  of  data  reduction  equations  is  found 
in  Appendix  D).  The  angles  are  taken  care  of  by  using  the  coordinate  transformation 
described  in  Appendix  B.  After  extreme  outliers  are  eliminated  (if  they  exist),  data  are 
reduced  once  to  determine  a,  then  again  with  user  specified  outlier  criterion  (±3a,  in  this 
study),  and  revised  statistics  are  calculated.  The  reduced  LDV  data  format  is  presented  in 
Appendix  E.  The  data  reduction  software  output  included  the  following:  coordinates  (x,  y,  z), 
reference  velocity,  normalized  mean  velocities,  normalized  turbulence  intensities, 
normalized  turbulent  stresses,  and  normalized  turbulent  triple  products.  Also,  the  raw  data 
from  the  LDV  signal  processor  are  stored  and  available  for  researchers  who  wish  to  obtain 
other  information. 
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CHAPTER  4 


EXPERIMENTAL  PROCEDURE 


4.1  Introduction 

The  mean  flowfield  parameters  behind  a  bluff  body  mounted  at  mid-stream  in  a 
rectangular  test  section  were  measured  using  both  a  1-d  and  3-d  laser  Doppler  velocimeter 
(LDV).  A  description  of  the  experimental  apparatus  used  to  make  these  measurements  was 
given  in  Chapter  3.  This  chapter  describes  the  techniques  used  to  obtain  the  following 
flowfield  parameters: 

,  1 .)  Mean  axial,  transverse,  and  spanwise  velocities  {U,V,W). 

2. )  Mean  axial,  transverse,  and  spanwise  rms  turbulence  levels  (m,  v,  w). 

3. )  Mean  spanwise  velocities  (Vf )  along  the  z-axis  at  y=0. 

4. )  Turbulent  normal  stresses  (mm,  vv,  ww). 

5. )  Reynolds  stresses  (mv, MW,  vw). 

'  6.)  Turbulent  triple  products  (mmm,  mmv,  mmw,  mvv,  uww,  wv,  ww,  vww,  www). 

Also  included  in  this  chapter  are  parameter  settings  for  the  various  LDV  systems  and 
flow  system  parameters,  as  well  as  the  statistical  uncertainty  estimation  techniques. 
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4.2  Test  Conditions 


As  noted  previously,  the  primary  objective  of  this  study  was  to  provide  an  accurate, 
experimentally  measured  data  set  for  the  turbulent  flow  parameters  behind  a  bluff  body 
mounted  in  a  rectangular  test  section.  With  this  in  mind,  all  flow  conditions  were  carefully 
monitored  and  documented  throughout  the  testing  procedure.  The  flow  conditions  also 
included  atmospheric  temperature  and  pressure.  The  total  temperature  was  sampled  from 
thermocouples  via  an  A/D  converter  board  inserted  in  a  personal  computer  to  within  ±0.5  °F. 
The  static  pressure  was  also  sampled  via  a  pressure  transducer  to  ±0.01  psia.  The  pressure 
differential  between  static  and  total  was  constantly  monitored  using  a  highly  sensitive 
manometer  and  was  measured  to  within  ±0.025  inch-H^O.  The  inlet  velocity,  was  then 

calculated  for  every  measurement  point  at  every  station,  with  an  uncertainty  of  ±0.26  m/s 
(see  Table  4.1  and  Appendix  G).  The  reference  velocity  changed  slightly  from  one  data  point 
to  another  and  from  one  day  to  another,  depending  upon  the  local  atmospheric  conditions. 
However,  all  parameters  were  normalized  point-to-point  by  the  measured  reference  velocity. 

Table  4.1.  Reading  Variables  Error. 


data 

uncertainty 

T. 

±0.5°  F 

P., 

±0.01  Psia 

AP 

±0.025  in-H,  O 

U.f 

dhO.26  m/s 
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4.3  Sampling  Method 


In  all  cases,  5120  individual  realizations  were  accumulated  for  each  velocity 
component  at  each  measurement  point  to  form  turbulence  statistics.  The  coincidence  timing 
logic  on  the  TSI  IFA-750  master  interface  ensured  that  simultaneous  (within  20  |is)  data 
were  measured.  As  mentioned  earlier,  the  velocity  measurements  presented  in  this  study 
were  not  corrected  for  velocity  bias.  Mass  flow  rates  based  on  integration  of  axial  velocity 
measurements  at  all  stations  in  the  far-field  were  found  to  be  within  ±2%  of  the  upstream 
profile  at  station  x/H=-5  (Figure  5).  These  axial  velocities  were  measured  from  top  to  bottom 
to  include  measurements  in  the  boundary  layers. 


-5  0  5  10  15  20 


Figure  5.  Normalized  mass  flux  variation  from  station  to  station  based  on  measured  axial 
velocities. 
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4.4  Turbulence  Statistics 


The  techniques  used  to  obtain  the  following  flowfield  parameters  are  explained 
below.  In  computing  the  statistical  parameters,  all  points  lying  far  outside  of  the  data  set,  are 
eliminated  first.  Then  any  point  lying  outside  ±3a  was  discarded.  On  the  average,  the 
number  of  points  that  were  eliminated  was  less  than  50.  Out  of  5120,  this  number  is 
considered  on  the  high  side,  typically  20  outliers  is  an  average  theoretical  number  for  a  two- 
component  LDV  system.  However,  this  is  not  an  unreasonable  number,  since  in  this  study, 
elimination  or  discrimination  is  based  on  three  components  of  velocities,  and  if  one 
component  is  eliminated,  the  two  other  corresponding  components  are  also  discarded.  After 
the  ±3a  cutoff  was  applied  to  the  5120  realizations,  the  revised  statistics  were  calculated. 
The  mean,  variance,  cross-product  and  triple-product  were  then  computed  using  the 
following  equations  (see  also  Appendix  D); 


m  =  j-Z{u,-uY 


xy=j^t(u,-utv,-v) 


uuu  = — ~V ) 
N  ,=i 


where  N  is  the  number  of  actual  samples  in  the  revised  data  set. 


(4.1) 

(4.2) 

(4.3) 

(4.4) 
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4.5  Statistical  Error  Estimation 

The  irregularity  of  the  turbulent  flowfield  necessitates  the  use  of  statistical  averaging 
methods  for  the  determination  of  turbulence  parameters.  Therefore  a  statistical  uncertainty 
^alysis  is  deemed  necessary  mainly  because  it  is  impossible  to  determine  exactly  any 
statistical  quantity  (random  variable)  using  a  finite  sample  size.  This  type  of  uncertainty  is 
based  on  estimation  and  sampling  theory.  For  instance,  the  standard  deviation  of  a  random 
distribution  of  mean  values  of  a  variable  F  is  expressed  as; 

(4.5) 

where  “0“  is  the  standard  deviation,  and  “N”  is  the  number  of  actual  samples  in  the  revised 
data  set. 

For  a  normal  distribution,  a  level  of  confidence  can  now  be  written  for  the  measured 
values  of  F  by  including  an  appropriate  coefficient.  For  instance,  for  a  95%  confidence 
interval  ((x/2=0.025),  the  equation  for  AF  is; 

AF  =  1.96a  (4.6) 

The  equation  states  that  the  computed  values  of  F  will  not  be  in  error  from  the  true 
ihean  value  by  more  than  ±AF,  95%  of  the  time. 

4.5.1  Uncertainty  in  the  Means 


The  uncertainty  equations  for  the  mean  velocities  are  expressed  as  follows; 


A(7  =  ±1.96a„/V}v 

(4.7) 

AV=±1.96a^/Viv 

(4.8) 

AW  =  ±1.96a^/Viv 

(4.9) 
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Using  these  equations,  for  a  95%  confidence  (Snyder,  et  al.  984),  the  maximum 
absolute  expected  error,  for  the  axial  velocity,  was  found  to  be  1.17%  at  station  x/H=0.5, 
where  the  axial  turbulence  intensity  level  was  42%.  The  maximum  absolute  expected  error 
for  the  transverse  velocity  was  found  to  be  1.92%  at  station  x/H=0.9,  where  the  turbulence 
intensity  level  was  70%.  As  for  the  spanwise  velocity,  the  maximum  absolute  expected  error 
was  found  to  be  0.95%  at  station  x/H=0.5,  where  the  turbulence  intensity  level  was  34%  (see 
Appendix  G). 

4.5.2  Uncertainty  in  the  Higher  Order  Moments 

In  this  study,  the  jackknife  resampling  method  described  by  Benedict  &  Gould  (1996) 
was  used  to  calculate  the  uncertainties  of  higher  order  statistics  such  as  the  2nd  and  3rd  order 
moments.  These  terms  included  the  six  turbulent  stresses,  and  all  nine  turbulent  triple 
products.  All  uncertainties  are  represented  by  error  bars  in  the  experimental  data  profiles. 
Also,  these  uncertainties  were  used  to  estimate  the  uncertainty  error  in  the  TKE  terms  (see 
section  6.9). 

4.6  Measurement  Errors 

When  using  LDV  as  a  tool  to  perform  experimental  measurements,  several  sources  of 
errors  might  result.  These  errors  usually  result  from  optical  misalignment,  such  as  the 
alignment  of  the  optical  table  and  the  test  section,  signal  processing,  refractive  index 
variations  of  the  gas  medium  and  of  the  windows,  probe  volume  effects,  particle  seeding 
effects,  and  the  data  sampling  method  used.  However,  most  of  these  errors  can  be  reduced  to 
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an  insignificant  level  or  even  eliminated.  For  instance,  in  this  study  the  angle  between  the 
test  section  and  the  optics  table  was  measured  and  found  to  be  less  than  0.03  degree. 

There  are  various  sources  of  biasing  error  that  can  effect  the  accuracy  of  LDV 
measurements  made  in  highly  turbulent  flows.  These  are  signal  bias,  fringe  bias,  and  velocity 
bias.  Signal  bias  (also  called  directional  bias)  occurs  mainly  in  recirculating  flows  such  as 
the  flow  studied  here.  Particles  traveling  in  the  opposite  direction  to  the  main  flow  cannot  be 
discriminated  from  particles  flowing  at  the  same  velocity  in  the  main  direction  by  the 
photodetector,  since  they  both  scatter  light  the  same  way.  However,  frequency  shifting 
imposes  a  fixed  frequency  on  the  fringes,  typically  40  MHz,  which  modifies  the 
photodetector  output  so  that  information  from  turbulent  and  reversing  flows  can  be  evaluated 
accurately.  Now,  the  fringes  move  so  that  a  zero  velocity  particle  will  produce  a  Doppler 
frequency  equal  to  the  shift  frequency.  Fringe  bias  occurs  when  seed  particles  traverse  the 
fringes  of  the  probe  volume  at  large  angles  of  incidence  and  generate  an  incomplete  number 
of  fringe  crossings  required  to  allow  signal  validation  by  the  processor.  This  reduces  the 
signal  validation  of  low  velocity  particles  which  leads  to  a  bias  towards  an  incorrectly  high 
value  of  mean  velocity.  However,  this  type  of  bias  error  can  also  be  eliminated  by 
appropriate  frequency  shifting  (Roesler  et  al.,  1980).  Velocity  bias  occurs  because  the 
velocity  data  are  obtained  as  the  seed  particles  randomly  pass  through  the  probe  volume. 
Naturally  this  passage  of  particles  takes  place  at  unequal  time  intervals.  Due  to  the 
fluctuating  nature  of  turbulent  flow,  there  will  be  more  particles  per  unit  time  passing 
through  the  probe  volume  during  time  intervals  when  the  velocity  is  high.  This  leads  to  a 
bias  towards  high  mean  velocity  as  calculated  by  simple  ensemble  averaging  of  the  data.  The 
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magnitude  of  this  error  frequently  outweighs  most  of  the  other  errors  listed  above  (Gould  et 
al.,  1986). 

The  techniques  currently  used  for  eliminating  velocity  bias  error  include  the  1-d  and 
2-d  corrections  of  McLaughlin  and  Tiederman  (1973),  and  the  time  between  validated  data 
method  of  Barnett  and  Bentley  (1974).  These  techniques  can  be  used  for  bias  correction  after 
the  data  are  acquired.  Velocity  bias  correction  can  also  be  done  during  data  acquisition. 
Experimental  sampling  techniques  have  also  been  proposed,  such  as  the  equal  time  interval 
sampling  technique  first  proposed  by  Stevenson,  et  al.  (1982).  However,  this  method 
requires  a  relatively  high  data  rate  on  the  order  of  three  times  the  passage  frequency  of  the 
turbulent  microscales  as  concluded  by  Gould,  et  al,  (1986).  It  was  found  that  the  highest 
passage  frequency  of  the  turbulent  microscales  in  a  turbulent  flow  behind  an  axisymmetric 
sudden  expansion  was  approximately  2000  Hz.  Therefore,  for  equal  interval  sampling  to 
work  for  such  a  geometry,  the  data  validation  rate  must  be  greater  than  6000  Hz.  After 
examination  of  post  facto  bias  correction  schemes,  Gould  et  al.  (1986)  concluded  that  these 
schemes  were  ineffective  and  could  introduce  larger  errors  than  the  original  uncorrected  data 
contained. 

The  coincident  data  rate  in  this  study  varied  from  5000  per  second  in  the  freestream, 
to  10(X)  per  second  behind  the  bluff  body,  in  the  recirculation  zone,  for  a  clean  window. 
Gould,  et  al.  (1986)  found  that  equal  interval  sampling  becomes  less  effective  as  the  data 
validation  rate  dropped,  to  the  point  where  no  bias  elimination  takes  place  as  the  data  rate 
approaches  500  Hz. 

In  this  study,  the  results  are  believed  to  be  slightly  biased  in  most  of  the  high 
turbulence  regions,  especially  in  the  recirculation  region.  However,  the  mass  balance 
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variation  from  station  to  station  shows  that  this  bias  error  in  the  regions  where  biasing  is 
possible,  accounts  for  a  global  mass  flow  error  of  less  than  2%.  Local  errors  in  the  mean 
velocities  are  probably  larger  than  this.  In  fact,  Gould,  et  al.  (1986)  gave  a  rule  of  thumb 
which  gives  the  error  in  mean  velocity  due  to  biasing  ~l/3  the  turbulence  intensity. 

Other  possible  error  sources  which  may  effect  an  LDV  measurement  by  reducing  the 
signal  to  noise  ratio,  are  those  cau.sed  by  dirty  windows,  non-flat  windows,  non- 
homogeneous  windows,  and  impurities  and  small  chips  in  the  window  material.  These  error 
sources  change  the  optical  path  length  which  can  deviate  the  incoming  laser  beams.  The 
beams  also  tend  to  diverge  and  lose  spatial  coherence  (Baker  et  al.,  1974),  causing  a  drift  in 
the  beam’s  coincidence,  and  reducing  the  effective  size  of  the  probe  volume.  Among  other 
causes,  the  signal-to-noise  ratio  decreases  as  a  result  of  the  improper  beam  crossing.  The 
high  quality  quartz  windows  (originally  used  for  Schlieren  measurement)  used  in  this  study 
and  the  continuous  cleaning  of  the  windows,  after  every  three  data  points  (in  the 
recirculation  region)  helped  to  maintain  the  data  rate  at  the  maximum  values  attainable. 
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CHAPTERS 


FAR-FIELD  EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 

5.1  Introduction 

In  order  to  map  the  flowfield  throughout  the  tunnel,  the  velocity  measurements  were 
made  in  two  separate  sets:  the  “far-field”  and  the  “near-field”.  The  near-field  experimental 
measurements  are  comprised  of  locations  between  x/H=0.25  and  1.2  and  are  presented  in 
Chapter  6.  The  far-field  measurements  included  stations  from  5  base  heights  (x/H=-5) 
upstream  to  22  base  heights  (x/H=22)  downstream  of  the  bluff  body.  The  bluff-body  dump 
plane,  i.e.,  the  trailing  edge,  was  assigned  station  x=0.  Negative  x-values  refer  to 
measurement  stations  upstream  of  the  bluff  body.  Care  was  taken  to  maintain  a  symmetric 
flowfield;  measurements  were  taken  across  the  entire  vertical  (i.e.  transverse  coordinate) 
tunnel  cross  section  to  verify  that  the  flow  was  symmetric.  All  measurements  are  normalized 
with  respect  to  the  reference  velocity  measured  at  the  same  instance,  and  the  bluff-body 
height,  H=12.7  mm  (0.5  in).  The  three-component  LDV  system  provided  the  axial, 
transverse,  and  spanwise  mean  and  fluctuating  components  of  the  velocity  field  together  with 
the  normal  and  Reynolds-shear  stresses,  and  the  higher  order  turbulent  triple  products.  The 
results  of  the  far-field  measurements  are  presented  next. 

Two  separate  sets  of  measurements  are  presented  in  this  section.  The  LDV  beam 
configuration  of  the  three-component  LDV  system  was  such  that  wall-to-wall  measurements 
of  the  three  velocity  components  was  not  possible.  For  this  reason,  it  became  necessary  to 
make  two  separate  sets  of  velocity  distribution  measurements.  First,  using  the  LDV  system 
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in  a  single  component  configuration,  wall-to-wall  (top  to  bottom)  measurements  of  the  axial 
velocity  distribution  (U)  were  made  at  stations  x/H=-5,  -3,  1,  1.5,  2,  3,  4,  6,  8,  10,  14,  18,  and 
22.  The  single  component  measurements  were  used  for  two  purposes:  1.)  to  validate  the 
three-component  measurements  and,  2.)  to  check  for  conservation  of  mass.  The  mass  flow 
rate  was  obtained  by  numerically  integrating  the  axial  velocity  profiles  at  various  stations 
across  the  tunnel.  The  mass  flow  rates,  shown  in  figure  5,  were  found  to  be  within  ±2%  of 
the  upstream  profile  obtained  at  xyH=-5.  Three-component  LDV  measurements  were  made  at 
the  same  axial  stations  listed  above.  However,  measurements  were  limited  to  the  central  5.58 
cm  (2.2  in)  portion  of  the  flowfield  in  the  vertical  (transverse)  direction  (-2.2<y/H<2.2). 

5.2  Mean  Velocity  Data 
5.2.1  1-D  LDV 

Figures  6  and  7  show  the  evolution  of  the  axial  mean  velocity  and  turbulence 
intensity  distributions  at  all  measured  stations,  as  measured  with  the  single  component  LDV. 
The  vertical  bars  in  the  following  figures  denote  the  statistical  uncertainty.  The 
measurements  show  a  relatively  flat  upstream  velocity  profile  representative  of  fully 
developed  channel  flow  at  x/H=-5.  Nearer  to  the  bluff  body,  at  x/H=-3  the  incoming  velocity 
distribution  was  affected  by  the  presence  of  the  flameholder  and  shows  a  slight  velocity 
deficit  in  the  center.  Downstream  of  the  bluff  body,  highly  turbulent  regions  clearly  identify 
the  wake  and  show  the  extent  of  mixing  and  the  flow  recovery.  However,  this  recovery  was 
very  slow  and  even  at  the  last  measurement  station,  x/H=22,  the  presence  of  the  bluff  body 
was  still  evident  in  the  wake. 
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Figure  6.  Normalized  mean  axial  velocity  profiles;  1  -d  LDV. 


x/H=1.5  x/H=2  >c/H=3  x/H=4  x/H=6  x/H=8  x/H=10  x/H=14  x/H=18  x/H=22 
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Figure  7.  Normalized  axial  turbulence  intensity  profiles;  1-d  LDV. 


Figure  8  shows  the  distribution  of  the  axial  velocities  and  the  turbulence  ittteilsities  at 
the  two  upstream  locations.  Upstream  of  the  bluff  body,  x/H=-5,  the  flow  exhibits  a  flat 
velocity  profile  typical  of  fully  developed  channel  flow  and  a  uniform  3%  turbulence 
intensity  distribution  in  the  center  region  increasing  to  a  maximum  of  8%  within  the  wind 
tunnel  wall  boundary  layer.  At  these  low  levels  of  turbulence  one  does  not  expect  to 
encounter  velocity  bias  (Gould  et  al.,  1986).  Measurements  at  x/H=-3  show  the  upstream 
influence  of  the  bluff  body  on  the  incoming  flow,  typical  of  subsonic  flow  conditions,  The 
central  velocity  deficit  results  from  the  stagnation  point  on  the  front  of  the  bluff  body. 
However,  no  change  was  noticed  in  the  turbulence  leVel.  The  measurements  at  x/H=-5  were 
used  as  the  initial  conditions  in  the  numerical  studies. 

5.2.2  3-D  LDV 

Axial  mean  velocity  and  turbulence  intensity  distributions  measured  at  x/H=l  with 
the  one  and  three-component  LDV  systems  are  compared  in  Figure  9.  This  is  the  most 
demanding  measurement  station  in  the  far-field,  where  the  turbulence  intensity  is  the  highest 
and  the  measurements  are  taken  within  the  unsteady  recirculating  wake  behind  the  bluff 
body.  The  general  agreement  between  the  two  measurement  systems  at  all  stations  was  also 
excellent.  However,  the  slight  discrepancy,  between  the  two  profiles  shown  in  Figure  9,  may 
be  due  to:  1.)  differences  in  atmospheric  relative  humidity  (see  Appendix  G)  between  the 
time  when  the  two  profiles  were  taken  (winter  vs.  summer);  and  2.)  the  change  in  the  bluff- 
body  geometry,  which  varied  slightly  during  the  course  of  the  research,  due  to  surface 
erosion  caused  by  the  seeding  material.  Three-component  velocity  measurements,  closer  to 
the  bluff  body  than  x/H=l,  are  presented  in  the  near-field  measurements,  in  Chapter  6. 
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Figure  8.  Normalized  mean  axial  velocity 
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Figure  9.  Normalized  mean  axial  velocity  and  turbulence  intensity  profiles;  1-d,  3-d  LDV. 


Measurements  at  x/H=l  showed  a  large  velocity  deficit  in  the  center  of  the  wake.  Because 
this  measurement  location  is  just  downstream  of  the  stagnation  point  (see  Figure  10),  no 
negative  mean  velocities  were  realized.  The  near-field  measurements  will  show  that  the 
recirculation  bubble  extended  only  to  x/H=0.90. 


Tunnel  Wall 
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Figure  10.  Bluff-body  near-field  flow  description. 


In  Figure  9,  the  velocity  distribution  revealed  two  symmetric  shear  layers,  each  with 
a  thickness  of  0.6H.  The  resulting  wake  at  this  station  was  20%  wider  than  the  bluff-body 
height.  The  maximum  turbulence  level  was  realized  within  the  shear  layers,  roughly  at 
y/H=±0.2,  where  the  mean  velocity  gradient  was  the  largest.  As  anticipated,  downstream  of 
x/H=±0.2,  turbulence  production  dropped  as  the  velocity  gradient  decreased  and  the  flow 
became  more  uniformly  distributed.  Although  the  measurements  closer  to  the  centerline,  -0.2 
<y/H<0.2,  showed  velocity  gradients  of  the  same  magnitude,  the  turbulence  intensity 
dropped.  Before  offering  an  explanation  for  this  behavior,  the  distribution  of  the  transverse 
velocity  and  turbulence  field  are  examined. 
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Transverse  mean  velocity  and  turbulence  intensity  profiles  for  station  x/H=l  are 
shown  in  Figure  11.  The  effect  of  the  bluff  body  becomes  readily  observable  in  this  figure, 
where  the  flow  rapidly  turns  from  both  sides  toward  the  centerline  to  fill  the  low  pressure 
region  behind  the  bluff  body.  The  turbulence  intensity  reaches  its  maximum  level, 
v/U^=0.75,  at  the  centerline,  y/H=0.0,  corresponding  to  the  location  with  the  highest  mean 
velocity  gradient.  Interestingly,  the  regions  of  maximum  turbulence  production  in  the  axial 
and  transverse  directions  do  not  correspond  spatially.  The  following  is  provided  as  a  possible 
explanation.  It  is  well  established  that  this  configuration  results  in  periodic  shedding  of 
vortices  from  the  upper  and  lower  separation  points.  Time-averaged  estimates  of  the  rms 
turbulence  fluctuations  obtained  with  an  LDV  system  result  from  the  superposition  of  small 
scale  turbulence  and  large  scale  fluctuations  due  to  the  passage  of  the  periodic  vortices 
through  the  probe  volume.  In  standard  configuration,  the  LDV  instrument  is  incapable  of 
differentiating  between  these  two  sources  of  turbulence.  Assuming  that  the  flowfield  is 
spatially  and  temporally  symmetric,  the  passage  of  the  upper  and  lower  vortices  skew  the 
axial  velocity  probability  density  function  (PDF)  in  a  symmetrically  similar  fashion  at  the 
centerline  y/H=0,  as  shown  by  the  histogram  in  Figure  12.  However,  away  from  the 
centerline,  this  contribution  is  dissimilar,  i.e.,  different  regions  of  the  upper  and  lower 
vortices  pass  through  the  probe  volume  and  result  in  compilation  of  skewed  velocity  PDFs 
and  this  longer  moment  arm  relative  to  the  mean.  Hence,  a  higher  axial  turbulence  intensity 
is  measured  away  from  the  center  as  shown  in  Figure  9.  In  contrast,  the  contributions  of 
periodic  vortex  shedding  from  the  upper  and  lower  separation  points  to  the  transverse 
velocity  PDFs  are  directional.  For  example,  the  time-averaged  value  of  the  transverse 
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Figure  1 1 .  Normalized  mean  transverse  V( 
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velocity  at  the  y/H=0  is  negligible,  i.e.,  V  =0.  Therefore,  the  passage  of  a  vortex  shed  from 
the  upper  corner  of  the  bluff  body  skews  the  velocity  PDF  to  show  a  negative  (lower)  mean 
velocity  and  conversely,  a  vortex  shed  from  the  lower  comer  skews  the  velocity  PDF  to 
show  a  positive  (higher)  mean  velocity.  Over  time,  the  combined  effect  results  in  a  much 
wider  bi-modal  velocity  PDF,  as  shown  in  Figure  13,  and  a  higher  estimated  value  of 
turbulence  intensity  regardless  of  the  measurement  location,  y/H.  The  trend  of  lower  axial 
turbulence  intensity  at  y/H=0  as  compared  to  the  levels  in  the  shear  layers  is  preserved  out  to 
station  x/H=3.  Downstream  of  x/H=3,  the  upper  and  lower  shear  layers  do  not  show 
individual  characteristics  and  have  merged  to  form  a  single  highly  turbulent  region  (Figure 
7).  The  need  for  mixing  augmenters  is  clearly  evident.  For  example,  27.94  cm  (11  inches) 
downstream  of  the  bluff  body,  at  x/H=22,  the  shear  layers  have  only  grown  to  4  times  the 
bluff-body  thickness  and  occupy  roughly  67%  of  the  tunnel  cross  sectional  area.  However, 
one  must  remember  that  this  study  is  taking  place  in  a  confined  flow. 

Due  to  the  two-dimensional  nature  of  the  bluff  body,  the  spanwise  velocity 
distributions  (W)  show  negligible  values  throughout  most  of  the  measurement  domain. 
However,  the  spanwise  velocity  distribution  at  x/H=l,  shown  in  Figure  14,  suggests  that  the 
flowfield  immediately  downstream  of  the  bluff  body  is  three-dimensional.  After 
investigating  the  flowfield  in  the  spanwise  direction,  it  was  deduced  that  the  small  positive 
w-component  velocity  value  at  y/H=0  could  be  due  to  a  boundary  layer  asymmetry.  It  is 
believed  that  the  combined  effects  of  side  wall  boundary  layers  and  the  pressure  field  behind 
the  bluff  body  drive  the  flow  away  from  the  walls  and  towards  the  center  of  the  flameholder 
resulting  in  the  formation  of  two  pairs  (upper  and  Ipwer)  of  counter  rotating  spanwise 
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Figure  1 4.  Normalized  mean  spanwise 


vortices.  The  peak  value  of  turbulence  intensity  is  believed  to  be  due  to  the  unsteady  lateral 
motion  of  the  vortex  pairs. 

Figure  15  shows  a  comparison  between  the  three  measured  normalized  turbulence 
intensities  at  x/H=l.  Peak  values  of  axial  normalized  turbulence  intensity  (i.e.  the  square  root 
of  the  turbulent  normal  stress,  normalized  by  the  reference  velocity)  occur  at  x/H=l,  in  the 
shear  layers  and  reach  a  value  of  31%  of  the  freestream  mean  velocity.  This  value  is  slightly 
above  that  obtained  behind  a  backward-facing  step  (25%),  as  reported  by  Gould,  et  al. 
(1990).  It  is  noteworthy  that  the  peak  transverse  normalized  turbulence  intensity  is  77%, 
approximately  2.5  times  the  value  of  the  peak  axial  normalized  turbulence  intensity  (33%), 
and  2.75  times  that  of  the  spanwise  turbulence  intensity  (28%).  Axial  mean  velocity  and 
turbulence  intensity  profiles  obtained  with  the  three-component  LDV  system  are  shown  in 
Figures  16  and  17,  respectively.  There  exists  a  very  close  match  to  the  measurements 
obtained  with  the  single  component  LDV  system. 

Distribution  of  the  transverse  mean  velocities  and  turbulence  intensity  profiles  at  all 
measured  stations  are  presented  in  Figures  18  and  19,  respectively.  Negative  values  of  y/H 
correspond  to  measurement  locations  below  the  bluff-body  centerline.  At  x/H=-5,  the  inlet 
velocity  profile  shows  negligible  transverse  velocity  magnitude,  however,  at  x/H=-3,  there  is 
a  distinct  transverse  velocity  distribution.  Note  that  negative  values  of  transverse  velocity 
show  that  the  flow  has  begun  to  turn  downward  and  conversely  the  positive  values  indicate 
an  upward  turning  of  the  flow  around  the  bluff  body.  The  flow  rapidly  turns  to  fill  the  low 
pressure  wake  region  immediately  downstream  of  the  bluff  body.  The  peak  values  of  the 
transverse  velocities  are  located  inside  the  recirculation  zone  bound  by  the  upper  and  lower 
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Figure  15.  Axial,  transverse  &  spanwise 
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Figure  16.  Normalized  mean  axial  velocity  profiles;  3-d  LDV. 
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Figure  17.  Normalized  axial  turbulence  intensity  profiles;  3-d  LDV. 
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Figure  1 8.  Normalized  mean  transverse  velocity  profiles;  3-d  LDV. 
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Figure  19.  Normalized  transverse  turbulence  intensity  profiles;  3-d  LDV. 


shear  layers.  Comparisons  of  Figures  16  and  18  show  that  flow  recovery  is  much  faster  in  the 
fransverse  direction  than  in  the  axial  direction  with  essentially  zero  mean  transverse  velt)city 
downstream  of  x/H=3.  In  contrast,  the  axial  velocity  distribution  shows  a  wake  region 
characterized  by  a  velocity  deficit  that  continues  to  station  x/H=22. 

The  measurements  showed  virtually  no  spanwise  mean  velocity  as  shown  in  Figure 
20,  except  for  a  small  peak  at  x/H=l.  However,  the  normalized  spanwise  turbulence 
intensities  shown  in  Figure  21  exhibit  a  maximum  of  30%  at  station  x/H=l.  These  turbulent 
fluctuations  were  still  noticeable  in  the  wake  and  remained  until  station  x/H=18.  More  about 
this  will  be  discussed  in  the  near-field  section. 

5.3  Vortex  Shedding  Frequency 

In  an  attempt  to  find  the  vortex  shedding  frequency,  several  single  component 
measurements  with  large  sample  populations  at  high  seeding  rates  were  taken.  Since  LDV 
data  do  not  provide  a  uniform  time  interval  between  velocity  realizations,  direct 
transformation  of  the  velocity  data  from  the  time  domain  to  the  frequency  domain  is  not 
easily  accomplished.  In  order  to  convert  the  non-uniformly  spaced  velocity  data  into  uniform 
time  intervals,  an  averaging  technique  was  employed.  Essentially,  a  time  interval.  At,  was 
selected  and  the  data  falling  within  each  time  interval  were  averaged.  This  averaging 
technique  (time  history  method;  see  TSI  manual,  1992)  was  applied  to  the  complete  data  set 
to  achieve  an  approximate  velocity  vs.  equidistant  time  between  samples.  A  Fast  Fourier 
Transform  (FFT)  routine  was  then  used  for  direct  computation  of  the  frequency  spectrum. 
Figure  22  shows  the  power  spectrum  associated  with  velocity  samples  collected  at  x/H=-5 
and  x/H=2  at  the  indicated  y-locations.  A  distinct  frequency  peak  is  noted  at  roughly  /=454 
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Figure  20.  Normalized  mean  spanwise  velocity  profiles;  3-d  LDV! 
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Figure  22.  Frequency  Spectra  at  x/H=-5  &  x/H=2. 


Hz  at  the  x/H=2  location.  Power  spectra  of  measurements  at  other  downstream  locations 
repeatedly  identified  this  frequency,  except  at  x/H=22,  which  did  not  show  a  dominant 
frequency.  This  shedding  frequency  corresponds  to  a  Strouhal  number,  St=0.24,  where; 

(5.1) 

Power  spectra  of  the  upstream  velocity  measurements  did  not  identify  a  high  energy  peak  at 
any  particular  frequency. 

5.4  Turbulent  Normal  and  Shear  Stresses 

Figures  23,  24  and  25  show  profiles  of  the  three  measured  normalized  turbulent 
normal  stresses  (mm,w,ww).  Once  again,  the  bars  in  these  figures  denote  95%  confidence 
intervals.  The  axial,  transverse  and  spanwise  turbulent  normal  stresses  are  clearly  different  in 
value  and  profile  shapes.  The  different  shapes  indicate  that  the  structure  in  and  near  the 
recirculation  zone  is  quite  anisotropic  and  remains  so  until  x/H=8.  The  normalized  transverse 
turbulent  normal  stresses  exhibit  the  highest  peak  values  (0.6),  approximately  five  times  the 
value  of  the  peak  axial  turbulent  normal  stresses  (0.12).  This  is  due  primarily  to  estimate 
average  velocities  which  include  both  large  scale  periodic  structures  and  the  smaller  scale 
more  random  like  turbulence.  The  measured  normalized  turbulent  shear  stresses  (uv,uw,vw) 
are  shown  in  Figures  26,  27  and  28.  These  measurements  indicate  that  the  Mv-shear  stress  is 
the  only  significant  one  of  the  three  with  peak  normalized  values  of  ±0.12,  occurring  at 
x/H=l. 
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Figure  23.  Normalized  axial  turbulent  normal  stresses. 
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Figure  24.  Normalized  transverse  turbulent  normal  stresses. 
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Figure  25.  Normalized  spanwise  turbulent  normal  stresses. 


x/H=1.5  x/H=2  x/H=3  x/H=4  x/H=6  x/H=8  x/H=10  x/H=14  x/H=18  x/H=22 


Normalized  uv  turbulent  shear  stresses 


x/H=1.5  x/H=2  x/H=3  x/H=4  x/H=6  x/H=8  x/H=10  x/H=14  x/H=18  x/H=22 


84 


Figure  27.  Normalized  Ow  turbulent  shear  stresses. 


x/H=3  x/H=4  x/H=6  x/H=8  x/H=10  x/H=14  x/H=18  x/H=22 


Figure  28.  Normalized  vw  turbulent  shear  stresses. 


5.5  Turbulent  Triple  Products 


The  turbulent  triple  products  determined  from  the  measurements  made  in  this  study 
are  uuu,  uuv,  uvv,  uuw,  vww,  uww,  vw,  www,  and  ww;  and  are  shown  in  their 
normalized  form  in  Figures  29  through  37.  The  turbulent  triple  products  represent  the 
transport  of  the  turbulence  (TKE)  by  the  fluctuating  velocities.  The  maximum  normalized 
(absolute)  values  in  their  order  of  magnitude  are  wv/f/^^=0.28,  Mw/f/^^^=0.05, 

uuv  I  Ulf=0.03,  uuu  /  Ul^=0.02S,  vww =0.0 18,  www /t/ ’^=0.01,  vvw/C/^^^  =0.005, 
uuwlUl^=OW3,  and  www /C/^^  =0.0028.  The  bi-modal  hump-like  feature  of  mmm -triple 
product  is  due  to  the  upper  and  lower  shear  layers.  The  value  of  vw  was  the  largest  («10 
times  that  of  uuu).  For  the  same  reason  described  earlier  when  comparing  uv  and  uu .  It  is 
worth  noting  that  all  terms  containing  an  odd  power  of  v  such  as  uv ,  vw ,  etc.  are  anti¬ 
symmetric,  while  those  containing  an  even  power  of  v  such  as  uvv,  etc.  are  symmetric  about 
the  centerline.  This  is  mainly  due  to  the  dominant  antisymmetric  nature  of  the  transverse 
velocity  component  (Figure  18). 
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Figure  29.  Normalized  uuu  triple  correlations. 
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Figure  30.  Normalized  vw  triple  correlations. 
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Figure  31 .  Normalized  www  triple  correlations. 


Figure  32.  Normalized  uw  triple  correlations. 
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Figure  33.  Normalized  vww  triple  correlations. 
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Figure  34.  Normalized  vvw  triple  correlations. 
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Figure  35.  Normalized  uww  triple  correlations. 
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Figure  36.  Normalized  uuw  triple  correlations. 
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Figure  37.  Normalized  uuv  triple  correlations. 


^  _  MW  +  VV  +  WW  2) 
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Figure  38  shows  the  normalized  turbulent  kinetic  energy  profiles  at  the  measured 
stations.  Appendix  H  shows  how  the  statistical  uncertainty  for  k  is  calculated.  The  TKE 
peaked  at  station  x/H=l  at  the  centerline  to  a  normalized  value  of  0.38.  The  three  normal 
stresses  normalized  with  k  are  shown  in  Figures  39,  40  and  41.  These  figures  clearly  show 
that  the  turbulence  downstream  of  the  bluff  body  is  not  isotropic  and  thus  care  must  be  taken 
when  one  is  attempting  to  numerically  simulate  these  types  of  flow  using  simple  turbulence 
models.  Vortex  shedding  from  the  upper  and  lower  separation  points  is  believed  to  be  the 
main  cause  of  this  phenomenon.  The  isotropic  normalized  values  for  uu,  vv  and  ww  should 
be  equal  to  2k/3.  However,  such  a  value  was  found  only  in  the  freestream  or  at  further 

downstream  positions  behind  the  bluff  body.  The  value  for  vv  /  k  was  found  to  be  the 
highest  with  a  maximum  value  of  1.7,  while  that  of  ww  !  k  was  the  lowest  with  a  maximum 


value  of  1.1. 

It  is  obvious  from  the  measurements  shown  in  the  above  figure  that  there  is  much 
more  to  be  investigated  in  the  very  near-field  just  downstream  of  the  bluff  body.  This  will  be 
the  subject  of  the  next  Chapter. 
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Figure  39.  Axial  mean  normal  stresses-  Isotropy  check. 


x/H=1.5  x/H=2  x/H=3  x/H=4  x/H=6  x/H=8  x/H=10  x/H=14  x/H=18  x/H=22 


99 


Figure  40.  Transverse  mean  normal  stresses-  Isotropy  check. 
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Figure  41 .  Spanwise  mean  normal  stresses-  Isotropy  check. 


CHAPTER  6 


NEAR-FTFJ.D  EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 

6.1  Introduction 

In  this  section,  the  results  of  laser  Doppler  velocimeter  (LDV)  velocity  measurements 
in  the  “near-field”  are  presented.  Measurements  of  the  axial,  transverse  and  spanwise  mean 
velocities,  turbulent  normal  stresses  as  well  as  Reynolds  stresses  and  turbulent  triple 
products  are  included  here.  The  turbulent  kinetic  energy  (TKE)  as  well  as  the  turbulence 
dissipation  rate  were  also  calculated.  In  this  part  of  the  study,  the  blue  and  green  laser  beams 
were  rotated  to  ±45°  relative  to  the  tunnel  bottom  wall  in  order  to  gain  access  to  the  closest 
possible  location,  near  the  trailing  edge  of  the  bluff  body,  without  beam  interference.  The 
following  locations  were  investigated  in  the  near-field,  within  and  near  the  recirculation 
zone;  x/H=0.5,  0.6,  0.7,  0.8,  0.9,  1.0,  1.1,  1.2.  The  smoothness  of  the  profiles  indicate  that 
the  data  rate  and  the  sample  size  (N=5120)  are  quite  sufficient  to  obtain  smooth  data,  even 
for  the  higher  order  moments.  The  mean  spanwise  velocities  along  the  z-axis  at  y/H=0  at 
various  x/H  stations  are  also  shown  in  this  Chapter.  Comparisons  between  commonly  used 
assumptions  for  ww,  uww,  and  vww  (when  third  component  data  is  not  available)  and  the 
experimental  measurements  of  these  quantities  are  made  and  presented.  Newly  derived 
expressions  for  these  turbulence  parameters  are  suggested  using  the  data. 
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6.2  Mean  Velocity  Data 


Simultaneous  three-component  measurements  of  mean  velocity  distributions  are 
obtained  by  traversing  the  flowfield  with  the  LDV  probe  volume  in  the  vertical  (y)  direction 
at  various  axial  positions.  The  measurements  were  made  only  out  to  y/H=±l.2  as  the 
velocities  and  higher  order  moments  were  nearly  uniform  in  the  freestream  away  from  the 
centerline  region,  as  depicted  by  the  previously  presented  far-field  measurements. 

Non-dimensional  mean  axial  velocity  profiles  are  shown  in  Figure  42.  All  velocity 
measurements  are  normalized  with  the  reference  velocity,  as  obtained  p6int-by-point,  and  all 
dimensions  are  normalized  with  the  bluff-body  height,  H=12.7  mm  (0.5  in.).  Measurements 
gave  virtually  identical  results  when  compared  to  those  obtained  with  the  previous  LDV 
setups  (i.e.,  see  1-d  and  3-d  measurements  with  0°  and  90°  beam  orientation  at  x/H=l).  At  all 
stations,  the  presence  of  the  wake  behind  the  bluff  body  can  clearly  be  seen  in  the  center  of 
the  flow,  as  well  as  the  presence  of  two  shear  layers.  Downstream  of  the  bluff  body,  a  time- 
averaged  recirculation  zone  is  present  with  associated  reverse  mean  velocities.  The  first 
station,  located  one-half  of  the  bluff-body  height  downstream  (x/H=0.5),  shows  a  maximum 
normalized  velocity  deficit  of  -0.35  (35%  of  the  reference  velocity).  Although,  the 
recirculation  zone  is  often  called  the  “dead  air  zone,”  in  this  study  it  contained  unsteady 
flow.  The  length  of  this  ellipsoidal-shaped,  time-averaged  recirculation  zone,  behind  the 
bluff  body,  defined  where  t7  =  0  at  y/H=0,  was  found  to  be  x/H=0.9.  The  width  of  the  time- 
averaged  recirculation  zone  was  nearly  one  bluff-body  thick,  which  is  relatively  small,  when 
compared  to  the  recirculation  region  behind  a  backward-facing  step.  These  dimensions  are  in 
agreement  with  Raghunathan  (1994).  The  shorter  length  of  the  recirculation  zone  is 
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Figure  42.  Normalized  mean  axial  velocity  profiles;  3-d  LDV. 


attributed  to  the  inverse  relationship  that  exists  between  blockage  and  the  size  of  the 
recirculation  zone.  Davies  and  Beer  (1970)  indicated  that  the  ratio  of  the  pressure  forces  to 
inertia  forces  in  the  wake  of  the  bluff  body  increases  with  the  blockage  ratio,  resulting  in  a 
shorter  distance  to  the  stagnation  point.  However,  the  magnitude  of  the  reverse  flow  is  high 
in  this  region  which  results  in  rather  intense  vortex  role-up.  The  width  of  the  unsteady  wake 
extended  from  -0.6<y/H<0.6  for  all  eight  near-field  stations,  which  is  slightly  larger  than  the 
bluff-body  thickness.  The  normalized  mean  axial  velocities  at  y/H=0  increased  in  the 
downstream  direction  from  a  value  of  -0.35  at  x/H=G.5  to  0.6  at  x/H=1.2  (Figure  42). 

Figure  43  shows  the  normalized  transverse  mean  velocity  distributions.  The  mean 
transverse  velocity  profiles  are  anti-symmetric  at  all  stations,  since  the  transverse  velocity  is 
directed  towards  the  centerline  in  this  near-field  region.  Unlike  the  axial  velocity,  the 
magnitudes  of  the  transverse  velocities  were  relatively  unchanged  from  axial  station  to 
station  due  to  transverse  direction  oscillating  shed  vortices  which  span  the  entire  length  of 
the  time-averaged  recirculation  zone.  This  is  a  result  of  the  vortex  roll-up  process,  which  acts 
both  upward  and  downward  during  the  sampling  time.  In  the  axial  direction,  the  vortex 
contribution  to  the  U-velocity  decreases  as  the  convective  freestream  flow  begins  to 
dominate  the  weakening  recirculating  flow  with  increasing  downstream  distance. 

Figure  44  shows  the  normalized  mean  spanwise  velocity  (W)  measurements  at  z/H=0 
in  the  transverse  direction.  The  data  show  that  there  exists  a  w-velocity  component  at  the 
centerline  station  which  diminishes  in  magnitude  as  the  flow  progresses  further  downstream, 
away  from  the  bluff  body.  In  a  perfectly  symmetric  flow,  one  would  not  expect  such 
behavior;  W  would  be  zero  in  this  case.  These  measurements  may  indicate  that  the  probe 
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Figure  44.  Normalized  mean  spanwise  velocity  profiles  at  z/H=0;  3-d  LDV. 


volume  was  not  centered  properly  with  respect  to  the  fluid  flow  centerline  in  the  spanwise 
direction.  It  is  also  speculated  that  non-symmetric  boundary  layer  growth  on  the  side  walls 
could  lead  to  a  slight  shift  in  the  flow  centerline  in  the  spanwise  direction. 

Figure  45  shows  the  normalized  mean  spanwise  velocity  measurements  along  the  z- 
axis  at  y/H=0  in  the  near-field.  This  figure  clearly  shows  the  three-dimensional  nature  of  the 
flow  behind  the  bluff  body.  However,  this  w-component  is  not  symmetric  with  respect  to  the 
z/H=0  axis  but  shows  a  positive  value  for  W  as  was  shown  in  Figure  44.  In  addition,  Figure 
45  suggests  that  there  exists  four  counter-rotating  (in  the  x-z  plane)  recirculating  bubbles 
directly  behind  the  bluff  body  that  were  slightly  non-symmetric  and  thus  did  not  cancel  at 
z=0.  This  non-symmetric  behavior  could  be  attributed  to  many  factors,  such  as,  imperfect 
walls,  a  slight  shift  in  bluff-body  geometry,  non-symmetric  boundary  layers  as  mentioned 
earlier,  and  possibly  other  factors.  Three-dimensional  flow  characteristics  were  also  reported 
by  Armaly,  et  al.  (1983)  in  the  recirculation  region  behind  a  backward-facing  step. 

In  an  effort  to  define  the  recirculation  zone  location  in  the  spanwise  direction, 
separate  1-d  LDV  measurements  were  made  at  fifteen  closely  spaced  axial  locations  along 
the  z-axis.  As  shown  in  Figure  46  the  incoming  spanwise  flow  profile  is  flat  as  depicted  by 
the  two  upstream  stations.  Also,  downstream  of  the  bluff  body,  a  similar  general  profile 
shape  as  that  obtained  when  using  the  3-d  LDV  (i.e..  Figure  45)  is  noticed.  One  difference 
however,  is  that  the  sign  of  the  IV -velocity  trend  is  reversed  in  the  two  very  near  stations; 
x/H=0.02,  and  x/H=0.1,  when  compared  with  the  trend  of  the  other  downstream  stations. 
This  indicates  that  a  vortex  exists  downstream  of  the  bluff-body  trailing  edge  with  its  center 
located  very  near  the  bluff  body. 
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Figure  45.  Normalized  mean  spanwise  velocity  profiles  at  y/H=0;  3-d  LDV. 
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Figure  46.  Normalized  "1-d  LDV"  mean  spanwise  velocity  profiles  at  y/H=0;  1-d  LDV. 


6.3  Turbulent  Stresses 


6.3.1  Reynolds  Stresses 

The  normalized  turbulent  shear  stresses  (uv,uw,vw)  have  also  been  obtained  in  the 
near-field  region.  These  measurements  indicate  that  the  uv  shear  stress  (Figure  47)  is  the 
most  significant  one  of  the  three,  and  its  maximum  value  is  five  times  larger  than  the 
maximum  value  of  the  other  two  shear  stresses  (Figures  48  &  49).  The  values  of  uv  were 
found  to  be  rather  small  at  x/H=0.5,  the  closest  measurement  station  downstream  of  the  bluff 
body.  In  the  separated  region,  the  shear  stress  increases  to  an  absolute  maximum  of  about 
OAxU,^/  occurring  at  x/H=0.9,  y/H=±0.2,  which  is  the  locus  of  the  time-averaged 
separating  streamlines  between  the  recirculation  bubble  and  the  surrounding  flow. 

6.3.2  Normal  Stresses 

Figures  50  through  52  show  the  normalized  turbulent  normal  stresses  at  the  eight 
near-field  axial  stations.  It  should  be  noted  that  local  peak  intensities  occur  at  locations 
where  the  mean-velocity  gradients  are  high,  indicating  that  gradient  transport  turbulence 
models  may  work  well.  As  expected,  the  axial  normal  stresses  (Figure  50)  show  the  double 
hump  behavior  due  to  the  presence  of  shear  layers  near  the  top  and  bottom  edge  of  the  bluff 
body.  Peak  values  of  normalized  axial  turbulent  stress  occur  at  x/H=0.5,  in  the  shear  layers, 
and  reached  a  value  of  0.18,  which  corresponds  to  a  maximum  normalized  axial  turbulence 
intensity  of  42.5%.  The  normalized  value  (31%),  at  x/H=l,  compares  well  with  that  obtained 
behind  a  backward-facing  step  (25%),  as  reported  by  Gould,  et  al.  (1990).  It  is  noteworthy 
that  the  peak  normalized  transverse  normal  stress  occurs  on  the  centerline  of  the  flow  at 
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Figure  47.  Normalized  uv  turbulent  stress  distributions. 
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Figure  48.  Normalized  Ow  turbulent  stress  distributions. 
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Figure  49.  Normalized  vw  turbulent  stress  distributions. 
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Figure  50.  Normalized  uu  turbulent  stress  distributions. 
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Figure  51 .  Normalized  w  turbulent  stress  distributions. 
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Figure  52.  Normalized  wW  turbulent  stress  distributions. 


x/H=0.9,  and  is  roughly  0.49  (Figure  51).  In  addition,  the  peak  value  of  the  normalized 
spanwise  normal  turbulent  stress  (Figure  52)  was  found  to  be  0.12;  at  station  x/H=0.5. 

Durao,  et  al.  (1988)  made  two-component  LDV  measurements  behind  a  square  cross- 
section  cylinder  in  a  water  tunnel.  A  comparison  of  their  results  reveals  that  their  flow  is 
raiore  isotropic  than  this  bluff-body  flow,  but  larger  values  of  turbulent  normal  stresses  are 
found  in  their  flow  as  compared  with  the  ones  presented  here.  This  may  be  due  to  the  sharp 
comers  at  the  leading  edge  of  the  square  cross-section  cylinder.  In  particular,  on  the 
centerline  at  x/H=l,  w  =  2.4mh  in  their  flow;  while,  vv  =  5.3mm  at  the  same  location  in  this 
flow.  The  maximum  value  of  normalized  axial  turbulent  stress  on  the  centerline,  at  x/H=l, 
was  0.22  in  Durao,  et  al.’s  study  and  0.075  in  this  flow.  Similarly,  the  maximum  value  of 
normalized  transverse  turbulent  stress  on  the  centerline,  at  x/H=l,  was  0.68  in  their  study  and 
0.56  in  this  study. 

With  these  comparisons  in  mind,  one  must  note  that  the  flow  behind  a  bluff  body  is 
unsteady  and  somewhat  periodic  as  vortex  pairs  are  shed  from  the  separation  points.  As 
mentioned  in  Chapter  5,  a  dominant  frequency  of  454  Hz  (St=0.24)  was  found  for  this  flow. 
In  addition,  the  bi-modal  PDF  presented  earlier  (Figure  13)  gives  further  proof  that  the  flow 
contains  large  periodic  structures.  This  vortex  passage  is  responsible  for  the  large  peak 
values  of  transverse  turbulent  normal  stress  as  simple  ensemble  averaging  was  used  to  obtain 
the  turbulence  statistics.  It  is  believed  that  if  Durao,  et  al.’s  velocity  decomposition  method, 
which  separates  the  large  scale  turbulence  from  the  random  turbulence,  was  applied  to  this 
data,  this  flow  would  appear  more  isotropic,  as  resulted  in  Durao,  et  al.’s  study.  It  is 
noteworthy  that  the  dominant  frequency  behind  the  square  cross-section  cylinders 
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corresponds  to  a  Strouhal  number  of  0.13.  Again,  this  may  be  due  to  the  upstream  corners, 
which  may  produce  a  larger  recirculation  zone  thickness  than  the  actual  case. 

As  shown  by  the  far-field  measurements,  the  axial,  transverse  and  spanwise  turbulent 
normal  stresses  in  the  near-field  are  also  clearly  different  in  value  and  profile  shapes.  Again, 
the  different  shapes  indicate  that  the  structure  of  turbulence  in  and  near  the  recirculation  zone 
is  quite  anisotropic,  at  least  for  the  case  of  ensemble  averaging.  This  anisotropy  implies  the 
existence  of  multiple  length  scales  in  the  turbulence  structure. 

6.4  Turbulent  Kinetic  Energy  (TKE) 

Figure  53  shows  the  measured  normalized  turbulent  kinetic  energy  (TKE)  profiles  at 
the  near-field  stations.  The  error  bars  denote  the  95%  statistical  uncertainty  limits  (see 
Appendix  H).  A  center  peak  is  noticed  at  the  station  x/H=0.9,  which  corresponds  to  the 
stagnation  location.  This  maximum  TKE  occurs  on  the  centerline  behind  the  bluff  body  and 
has  a  normalized  value  of  32  percent.  Downstream  of  x/H=0.9,  this  single  peak  decreases  in 
magnitude,  while  two  off-centerline  peaks  are  depicted  in  the  recirculation  zone  upstream  of 
x/H=0.8.  This  center  deficit  region  has  a  maximum  normalized  value  of  16%  and  is  due 
mainly  to  the  influence  of  the  axial  normal  stresses  that  possess  a  double  hump  behavior  in 
the  recirculation  zone. 

6.5  Isotropy 

An  isotropic  flow  is  defined  as  one  with  no  average  shear  stress  and  therefore,  all 
three  normal  stresses  must  be  equal.  A  check  for  isotropy  was  performed  by  normalizing  the 
axial,  transverse,  and  spanwise  normal  stresses  with  the  turbulent  kinetic  energy,  and  are 
shown  in  Figures  54  through  56.  As  was  mentioned  in  the  previous  chapter,  isotropic 
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Figure  53.  Normalized  turbulent  kinetic  energy. 
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Figure  54.  Axial  mean  normal  stresses-  Isotropy  check. 
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Figure  55.  T ransverse  mean  normal  stresses-  Isotropy  check. 
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Figure  56.  Spanwise  mean  normal  stresses-  Isotropy  check. 


behavior  is  found  only  in  the  freestream  region  and  far  downstream  of  the  wake  region  in 
this  flow.  In  the  near-field,  anisotropy  is  clearly  evident  by  the  ratios  of  uu!  k  which  varied 
between  0.2  at  the  centerline  to  1 .0  in  the  shear  layers.  The  ratios  of  w  I  k  showed  a  peak 
value  of  1 .6  on  the  centerline,  while  ww!  k  showed  a  relatively  flat  value  of  0.3. 

6.6  WW  Approximation 

When  the  third  component  velocity  measurements  are  not  available,  researchers  rely 
on  assumptions  to  approximate  the  values  of  the  missing  terms.  For  hw  ,  the  following 
relationship  is  usually  used,  ww=C^{m+w)  where  C,=0.5  is  used  in  a  backward-facing  step 
geometry  (see  Driver  &  Seegmiller,  1985;  Fuji  &  Gomi,  1976).  The  turbulent  kinetic  energy 
(TKE)  defined  as  k  =  {uu  +  vv  +  wwy2  is  then  approximated  by  it=Q(MM+w)  where 

Cj=0.75  (see  Chandrsuda  &  Bradshaw,  1981;  and  Etheridge  &  Kemp,  1978).  It  should  be 
noted  that  for  the  case  of  isotropic  turbulence,  this  approximation  is  exact.  In  the  present 
study,  the  spanwise  velocity  terms  were  obtained  from  direct  measurements.  Validation  of 
the  estimation  of  ww  was  performed  by  direct  comparison  to  the  measured  data.  Figure  57 
shows  the  measured  and  estimated  normalized  turbulent  kinetic  energy  using  this  simple 
approximation.  The  dotted  lines  show  that  the  isotropic  approximation  clearly  overestimates 
the  actual  value  of  ww,  especially  in  the  recirculation  zone.  An  approximation  constant  of 
C,=0.2  (or  C2=0.6)  was  found  to  fit  the  data  of  this  study  much  better  than  the  value  derived 
from  isotropic  assumptions.  This  is  shown  in  Figure  57  by  the  continuous  solid  lines.  It  is 
noteworthy  that  a  single  constant  works  well  even  in  the  high  velocity  gradient  regions  of  the 
shear  layers.  Azad,  et  al.  (1987)  suggested  other  values  for  C,  which  varied  from  0.35  to  0.41 
for  fully  developed  pipe  flow,  from  0.33  to  0.48  for  flow  in  a  diffuser,  and  0.4  for  boundary 
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Figure  57.  Normalized  Turbulent  Kinetic  Energy  (Odata  k  approximated  by:  - - o.6(uu+w)  . o.75(mi+w)). 


layer  flow.  Values  for  C,  of  0.37  and  0.63  are  inferred  in  Table  5.2  from  Townsend  (1980) 
for  a  laboratory  equilibrium  boundary  layer  and  an  atmospheric  equilibrium  boundary  layer, 
respectively.  Andreopoulos  and  Wood  (1982)  found  that  a  constant  of  C,=0.5  did  not  match 
their  measured  constant  value  (as  reported  in  Figure  9)  in  the  turbulent  boundary  layer  with 
surface  roughness.  The  reason  that  a  small  value  for  C,  (=0.2)  better  matched  the 
experimental  data  of  this  study  was  the  high  anisotropy  encountered  in  this  flow.  In  order  to 
offset  the  high  values  of  vv,  which  were  approximately  a  factor  of  2.5  greater  than  the  uu 
values,  a  smaller  coefficient  is  needed. 

6.7  Turbulent  Triple  Products 

The  triple  velocity  correlations,  u^u-Ui^,  represent  the  diffusion  of  turbulent  stresses 

due  to  turbulent  velocity  fluctuations  and  are  called  the  turbulent  diffusion  flux  of  Reynolds 
stresses.  Their  gradients  constitute  the  diffusion  terms  of  the  turbulent  kinetic  energy  balance 
equation.  To  the  author’s  knowledge,  no  direct  measurement  of  these  terms  has  ever  been 
made  in  bluff-body  flows.  Researchers  have  proposed  various  models  to  simulate  these 
terms.  A  good  summary  of  current  models  is  given  by  Cormack,  et  al.  (1978)  where  a 
comparison  of  their  model  with  four  other  techniques  is  presented.  In  the  present  study,  all 
nine  turbulent  triple  products  {uuu,  uvv,  uww,  vuu,  vvv,  vvw,  wuu,  vww,  and  www)  were 
determined  and  are  presented  as  their  normalized  values  in  Figures  58  through  66.  The  error 
bars  denote  the  95%  statistical  uncertainty  limits.  These  uncertainties  are  proportional  to 
This  effect  can  be  seen  by  comparing  the  uncertainties  at  x/h=0.5,  y/h=-0.1  (Figures 
58-66)  and  other  neighboring  locations.  The  sample  size  at  this  one  location  was  erroneously 
set  to  1024  samples  instead  of  5120  samples.  The  profiles  of  uuu  had  a  double  hump  feature 
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Figure  58.  Normalized  uuv  turbulent  triple  product  distributions. 
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Figure  59.  Normalized  uuw  turbulent  triple  product  distributions. 
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Figure  61 .  Normalized  ww  turbulent  triple  product  distributions. 


x/H=0.5  x/H=0.6  x/H=0.7  x/H=0.8  x/H=0.9  x/H=1.0  x/H=1.1  x/H=1.2 


130 


Figure  62.  Normalized  uww  turbulent  triple  product  distributions. 
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Figure  66.  Normalized  www  turbulent  triple  product  distributions. 


with  negative  minimum  values,  in  contrast  to  the  uu  profiles.  This  double  hump  feature  of 
uuu  is  the  result  of  the  top  and  bottom  shear  layers  as  was  the  case  with  the  axial  turbulent 
normal  stress  distributions  shown  in  Figure  50.  The  profiles  of  vvv  were  antisymmetric, 
similar  to  the  V-profiles,  but  inverted.  The  mw  profiles  were  wavey  with  large  relative 
uncertainty,  and  no  particular  shape  could  be  detected.  It  is  noteworthy  that  the  absolute 
uncertainty  is  approximately  5  times  larger  for  vvv  (0.016  U]^)  than  for  www  (0.0035 

The  maximum  normalized  values  of  uuu  / and  wv  /  had  factors  of  5  and  3,  less  than 
the  normalized  values  of  uu/Ul^  and  respectively.  The  maximum,  value  of  vvv 

was  found  to  be  the  largest  of  all  the  turbulent  triple  products  (=5  times  that  of  uuu)  and 
peaks  at  the  shear  layer  locations  (y/H=±0.2)  and  at  the  stagnation  point  (x/H=0.9)  to  a 
normalized  maximum  absolute  value  of  0.2  (Figure  65).  All  transverse  component 
containing  triple  products  terms  showed  rather  small  values  in  the  recirculation  zone  and 
increased  to  maximum  values  at  the  stagnation  point;  otherwise  the  opposite  was  true. 

Once  again,  maximum  values  of  turbulent  triple  products  in  this  flow  were  found  to 
be  up  to  30  times  larger  than  those  reported  for  axisymmetric  sudden  expansion  (Gould  et 
al.,  1990)  and  the  backward-facing  step  flows  (Driver  and  Seegmiller,  1985).  This  implies 
that  the  ensemble  averaged  diffusion  of  turbulent  kinetic  energy  (which  includes  the  periodic 
vortex  shedding)  is  much  higher  in  this  flow  than  in  the  axisymmetric  sudden  expansion  and 
backward-facing  step  flows.  It  is  commonly  known  (Chandrsuda  and  Bradshaw,  1981)  that 
large  eddies  are  responsible  for  the  significance  of  the  triple  products.  The  maximum 
normalized  absolute  values  of  the  nine  triple  products  in  order  of  decreasing  magnitude  are: 
vvv =0.18,  uw  11)1^=0.06,  uuu!  Ul^=0.05,  =0.048,  vww/ 17^^^=0.016, 
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wwlUl^=O.Q\,  =0.008,  MMw/C/’^  =0.004,  =0.003.  Again  this  indicates 

that  the  iiw  triple  correlation  term  dominates  this  type  of  flowfield,  and  that  mixing  (or 
turbulent  diffusion)  is  highest  in  the  transverse  direction,  a  condition  highly  desirable 
because  it  allows  for  shorter  length  combustors. 

6.8  UWW  and  VWW  Approximations 

When  the  spanwise  velocity  component  cannot  be  measured,  the  following 
assumptions  are  generally  made  about  the  turbulent  triple  products  uww  and  vww : 
mvvw  =  C3(mw+mmm)  and  vhw=C4(mmv+vw),  where  C3=C4=0.50  (see  Driver  &  Seegmiller, 
(1985)).  Wood  and  Bradshaw  (1982)  obtained  a  value  of  C4=0.20  across  an  entire  turbulent 
mixing  layer  constrained  by  a  solid  surface.  Lawn  (1971)  deduced  from  Laufer’s  (1954)  fully 
developed  pipe  flow  results  that  vww  =  05vvv  =  05uuv  and  in  turn,  used 
r^  =  0.25(MMV  +  rw).  Nakagawa,  et  al.  (1975)  used  the  approximation,  vww=  vvv.  Azad,. 
et  al.  (1987)  deduced  that  values  for  ww  and  vww  vary  with  location  in  turbulent  flows  and 
that  generalization  is  very  risky.  These  approximations  for  uww  and  vww  may  have  a  large 
effect  on  the  turbulent  transport  of  the  turbulent  kinetic  energy.  Since  all  nine  turbulent  triple 
products  have  been  calculated  in  this  flow,  it  is  of  interest  to  determine  how  well  these 
suggested  approximations  predict  the  actual  values.  Figure  67  shows  a  comparison  between 
the  experimentally  determined  values  of  uww  and  some  of  the  approximations  mentioned 
above.  As  shown  by  the  dotted  lines,  a  value  of  €3=0.50  clearly  overpredicts  the  measured 
values.  There  is  a  great  mismatch  between  the  experimentally  determined  values  of  uww 
and  the  suggestions  of  Lawn  (1971)  and  Nakagawa  (1975).  No  single  value  for  C,  could  be 
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Figure  68.  Third  moment  vww  of  fluctuatuing  velocity.  (Odata  . -03=0.50 


found  for  this  flow,  but  as  shown  by  the  dashed  lines,  the  simple  relationship  uww  =0.2 uvv 
was  found  to  best  model  the  measured  values.  Figure  68  shows  a  comparison  between  the 
experimentally  measured  values  of  vww  and  the  commonly  used  approximations.  Again,  the 
modeled  approximations  greatly  overpredict  the  experimental  data  for  vww,  and  a  best  fit  to 
the  measured  values,  as  indicated  by  the  dashed  lines  in  Figure  68,  would  be  vww  =0.5mmv  . 
No  reasonable  value  for  could  be  defined  for  this  flow. 

6.9  Turbulent  Kinetic  Energy  fTKE)  Balance 


The  approximate  equation  for  the  transport  of  TKE  for  steady,  incompressible  flow  in 
Cartesian  coordinates  is  given  by  Equation  (6.1),  (also,  see  Appendix  F): 


dk  -dk 

U  —  +V  —  +  W  —  = 
ax  dy  dz 


(convection) 


1  d  uuu  d  uvv  d  uww  d  vuu  d  vvv  d  vww  d  wuu  d  ww  d  www 

2  dx  dx  dx  dy  dy  ^  dy  ^  dz  ^  dz  ^  dz 


(diffusion) 


—du  —dv  — dw  —du  —dv  —du  —aw  — av  —aw 

UU—-+VV—+WW——  +  UV——+UV——+UW——+UW +  vw —  +  vw 

dx  dy  dz  dy  dx  dz  dx  dz  dy 


(production) 

+£  (dissipation) 

Where  the  mean  axial,  transverse  and  spanwise  velocities  are  denoted  by  U,  V,  and  W 
respectively.  Equation  (6.1)  is  composed  of  the  following  terms:  1)  convection,  2)  turbulent 
diffusion,  3)  production,  and  4)  viscous  dissipation  of  turbulent  kinetic  energy;  where 
k  =  +  vv  +  h'iv)/2  .  The  exact  transport  equation  for  TKE  has  additional  unknown  terms 

which  are  not  included  in  Equation  (6.1).  These  include  the  three  pressure-diffusion  terms 
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which  were  assumed  to  be  negligible  in  this  study  (see  Appendix  F),  a  commonly  accepted 
assumption  at  locations  away  from  solid  walls.  The  viscous  dissipation  terms  are  represented 
by  £  (where  £  is  a  negative  value),  the  turbulence  dissipation  rate  (see  Appendix  F).  Since 
most  of  the  terms  given  in  Equation  (6.1)  were  measured  except  for  £,  this  equation  can  be 
used  to  estimate  the  turbulence  dissipation  rate  by  balance.  Therefore,  these  values  of 
dissipation  rely  heavily  on  the  accuracy  of  the  other  terms  in  Equation  (6. 1 ). 

Each  term  in  Equation  (6.1)  was  estimated  by  numerically  differentiating  the  data 
using  second  order  accurate  finite  difference  algorithms  (central  differencing  for  interior 
points  and  forward  differencing  for  end  points).  Since  it  is  well  known  that  numerical 
differentiation  is  a  noisy  process,  these  first  derivatives  of  the  data  were  then  smoothed  or 
least  squares  curve  fitted  to  produce  smooth  profiles  of  each  term.  The  smoothing  technique 
developed  by  Gould  (1994)  was  used.  The  selection  of  either  an  n*"  order  least  squares 
polynomial  curve  fit  or  the  Savitzky  and  Golay  (1964)  smoothing  algorithm,  which  applies  a 
convolution  kernel  to  the  original  data,  depended  on  the  shape  of  each  differentiated  data 
profile.  The  choice  of  smoothing  technique  was  based  on  how  well  each  worked;  the  best 
method  was  selected  for  each  individual  profile  of  each  variable  in  Equation  (6.1). 

Finally,  regions  in  the  flow  where  turbulent  transport  mechanisms  of  TKE  are 
dominant,  are  identified.  Figure  69a  shows  the  profiles  of  the  TKE  terms  at  the  bubble 
station,  x/H=0.9.  The  individual  terms  are  presented,  with  their  associated  statistical 
uncertainties  (see  Appendix  H)  in  Figures  69b  through  69e,  at  this  most  cntical  station. 
Statistical  uncertainty  calculations  at  all  other  stations  were  also  made,  but  all  uncertainties 
were  smaller  than  those  at  the  presented  station,  and  thus  are  not  shown.  Because  of  the  large 
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Figure  69a.  TKE  balance  at  station  x/H=0.9  (Bubble  Location). 
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Figure  69b.  Convection  term  with  error  bars  at  station  x/H=0.9  (Bubbie  Locatiorr): 
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Figure  69c.  Production  term  with  error  bars  at  station  x/H=0.9  (Bubble  Location). 
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Figure  69cl.  Diffusion  term  with  error  bars  at  station  x/H=0.9  (Bubble  Location). 
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Figure  69e.  Dissipation  term  with  error  bars  at  station  x/H=0.9  (Bubble  Location). 


errors  in  the  triple  products  caused  by  the  relatively  low  number  oft  samplfes  (5120),  the 
diffusion  term  shows  the  highest  uncertainty,  which  in  turn  gives  high  uncertainty  of  the 
dissipation  term,  since  this  last  term  is  obtained  by  balance. 

6.9.1  TKE  Terms 

Figure  70  shows  the  measured  values  of  the  total  convection,  production,  diffusion, 
and  dissipation  of  TKE  as  defined  by  Equation  (6.1)  at  various  axial  stations.  Statistical  ? 
uncertainty  bars  are  not  included  in  these  figures  for  clarity  sake  as  mentioned  earlier;  Also, 
note  that  “+”  denotes  Gain  and  denotes  Loss.  These  profiles  show  that  the  TKE  is  first 
produced  in  the  shear  layer  where  both  m  and  the  axial  mean  velocity  gradient,  dUldy,  are 
large.  This  gives  the  double  hump  profile  shown  in  this  Figure  at  x/H=0.5.  At  downstream 
locations,  however,  the  wdVfdy  and  uudU /dx  terms  contribute  to  the  production  of  TKE 
and  cause  the  profile  shape  to  change  from  the  double  hump  shape  to  the  single  peak  shape 
shown  in  Figure  69a  and  Figure  70.  Turbulent  diffusion,  which  is  difficult  to  model,  shows  a 
double  hump  feature  with  two  small  positive  peaks  at  the  shear  layers,  and  a  large  deficit  in 
the  center.  This  deficit  reaches  a  maximum  normalized  value  of  -0.062  at  the  stagnation 
point  (Figure  69a).  Convection  of  TKE  shows  a  double  hump  behavior  which  is  positive 
with  peaks  towards  the  outside  of  the  shear  layers,  and  is  near  zero  in  the  recirculation  zone 
behind  the  bluff  body  near  the  centerline.  The  convection  profiles  seem  to  maintain  constant 
values  and  shape  throughout  the  whole  near-field  region.  In  general,  the  data  show  that 
production  and  diffusion  significantly  contribute  to  the  turbulent  kinetic  energy,  with 
absolute  peaks  at  y/H=0.  Convection  remained  relatively  low,  while  diffusion  reached  a 
minimum  at  the  stagnation  point  (Figure  69a),  where  the  triple  products  are  high  indicating 
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Figure  70.  TKE  Balance. 
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Figure  71 .  TKE  Balance. 


the  presence  of  large  eddies  or  at  least  very  periodic  flow.  As  for  the  dissipation,  the  values 
remained  relatively  constant  and  were  of  the  same  size  as  the  diffusion  at  the  first  three 
stations.  Finally,  the  levels  of  convection,  diffusion,  production  and  dissipation  of  TKE 
diminished  relatively  quickly  in  this  flow  and  are  near  zero  for  x/H  >  3  as  shown  in  Figure 
71.  This  may  explain  why  the  profiles  of  TKE  (Figure  53)  relax  slower  in  the  axial  direction 
than  does  the  velocity  field;  TKE  is  not  diffusing  very  rapidly  at  these  downstream  locations. 

Figures  72  through  74  show  the  contributing  terms  to  each  of  the  production, 
convection,  and  diffusion  terms.  For  instance,  wdVfdy  was  the  dominating  term  in  the 
production  of  TKE  (Figure  72)  and  peaked  to  a  maximum  normalized  value  of  1,  at  the 
stagnation  station,  x/H=0.9.  The  Vdk  /  dy  term  was  the  only  one  that  had  a  significant  value 
in  the  convection  of  TKE,  with  a  maximum  double  peaks  of  normalized  value  of  0.25  at  the 
shear  layers.  In  the  diffusion  of  TKE,  the  term  d  wv/dy  shows  a  double  hump  behavior  with 
maximum  absolute  normalized  values  of  0.6  at  the  stagnation  station. 

6.9.2  Derivatives 

In  addition,  the  raw  measurements  were  numerically  differentiated  using  second 
order  accurate  finite  difference  algorithms  (central  differencing  for  interior  points  and 
forward  differencing  for  end  points)  and  plotted  in  Figures  75  through  81.  Note  that  no 
smoothing  was  performed  to  create  these  figures.  The  relative  smoothness  shown, 
particularly  for  the  1st  derivative  profiles,  give  additional  confidence  in  the  accuracy  of  the 
test  condition  and  the  measurement  data.  Additional  points  of  measured  data  (i.e.,  U,  uu, 
tiuu, ...)  and  first  derivatives  of  these  parameters  were  interpolated  so  that  smaller  Ax,,  could 
be  used,  giving  smaller  numerical  errors. 
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Figure  72.  Production  of  TKE. 
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Figure  73.  Convection  of  TKE. 
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The  X-  and  y-  derivatives  of  the  mean  axial  velocity  followed  the  same  decreasing 
trend  as  the  velocities  themselves  (Figure  75).  The  y-derivatives  showed  positive  and 
negative  peaks  at  the  shear  layers  while  the  x-derivatives  showed  a  maximum  peak  at  the 
centerline,  with  the  exception  of  the  two  first  stations  where  a  slight  dip  was  noticed.  The 
reverse  is  observed  in  the  case  of  the  transverse  derivatives  (Figure  76).  It  is  noteworthy  that 
the  x-transverse  derivatives  decreased  from  station  x/H=0.5  to  x/H=0.9  and  reversed  slope 
direction  as  they  began  to  increase  again  (Figure  76).  The  y-spanwise  derivatives  (Figure  77) 
showed  an  antisymmetric  behavior,  while  the  x-derivatives  reversed  direction  with  respect  to 
the  velocity  profiles. 

The  y-derivatives  of  the  mv -shear  stress  (Figure  78)  showed  minimum  centerline 
values  which  reached  a  maximum  at  the  stagnation  point.  The  first  station  was  the  only  one 
with  a  double  hump.  The  x-derivatives  showed  a  minimum  and  a  maximum  at  the  shear 
layers,  which  decreased  to  almost  zero  at  the  stagnation  point  and  began  to  increase  again. 
The  MM -derivatives  (Figure  79)  showed  a  decreasing  trend  away  from  the  bluff  body.  The 
transverse  y-derivatives  showed  a  maximum  and  a  minimum,  while  the  x-derivatives  showed 
a  double  hump.  The  transverse  normal  stress  (vv)  x-derivatives  (Figure  80)  showed  a 
minimum  and  a  maximum  which  decreased  away  from  the  bluff  body,  while  the  x- 
derivatives  showed  a  centerline  maximum  which  flipped  at  the  stagnation  point  and 
thereafter.  The  spanwise  normal  stress  x-derivatives  (Figure  81)  showed  an  anti-symmetric 
behavior,  while  the  y-derivatives  showed  a  reverse  direction  with  respect  to  the  velocity 
profiles. 
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Figure  75.  Normalized  mean  axial  velocity  and  velocity  gradient  distribution.s. 

(O  data,  steady  state  k-e, - unsteady  k-e, - unsteady  RNG  ) 
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Figure  76.  Normalized  mean  transverse  velocity  and  velocity  gradient  distributions. 

(O  data,  •  •  steady  state  k-e, - unsteady  k-e,  - unsteady  RNG ) 
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Figure  77.  Normalized  mean  span  wise  velocity  and  velocity  gradient  distributions. 
(O  data, - smoothed  curve-fit ) 
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Figure  78.  Normalized  uv  turbulent  stress  and  stress  gradient  distributions. 


Figure  79.  Normalized  uu  turbulent  stress  and  stress  gradient  distributions. 
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(O  data,  •  •  •  •  steady  state  k-£,  —  unsteady  k-e, - unsteady  RNG ) 
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Figure  81.  Normalized  ww  turbulent  stress  and  stress  gradient  distributions. 
(Odata,  ~ ^  smoothed  curve-fit) 


For  a  truly  two-dimensional  flow,  3U  ldx  =  -dV  I  By  (Etheridge  &  Kemp,  1977); 
however,  as  Figure  82  indicates,  this  flow  shows  a  departure  from  two-dimensionality,  that 
extends  even  beyond  station  x/H=1.2,  downstream  of  the  stagnation  point. 
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CHAPTER  7 


NUMERICAL  RESULTS  AND  DISCUSSION 

7.1  Introduction 

Like  the  experimental  studies,  many  CFD  researchers  have  conducted  numerical 
studies  on  simple  bluff-body  geometries.  These  studies  have  produced  solutions  using 
techniques  varying  in  sophistication  from  the  steady-state  algebraic  stress  model  of  Kodali 
and  Amano  (1986)  to  the  dynamic  subgrid-scale  Large-Eddy  Simulations  (LES)  of  Yang  and 
Ferziger  (1993).  Some  early  solutions  modeled  only  half  of  the  symmetric  geometry,  thus 
imposing  symmetry  on  the  flow.  However,  detailed  studies  such  as  the  current  study,  have 
shown  that  the  instantaneous  flowfield  consists  of  a  structured  pattern  of  alternating  large 
scale  vortices.  A  symmetry  plane  enforced  by  the  numerical  model  precludes  the  formation 
of  the  shed  vortices,  and  thus  requires  the  turbulence  model  to  account  for  this  physical 
mechanism.  Comparisons  to  experimental  measurements  have  shown  such  predictions  to  be 
inadequate  for  bluff-body  flowfields.  Researchers  have  been  able  to  improve  the  numerical 
modeling  by  including  the  entire  flowfield.  When  solved  with  time-accurate  algorithms,  the 
flow  is  free  to  oscillate  about  the  geometric  symmetry  plane.  Such  oscillations  add 
significantly  to  both  the  time-averaged  turbulent  normal  and  shear  stresses  in  the  wake  of  the 
bluff  body. 

7.2  Numerical  Technique 

The  numerical  analysis  for  the  current  study  was  executed  by  CFD  Research 
Corporation  (CFDRC)  of  Huntsville,  Alabama.  A  full  report  of  the  numerical  technique  and 
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results  can  be  found  in  Raffoul,  et  al.  (1995).  The  numerical  analysis  was  accomplished 
using  a  commercially  available  CFD  code,  CFD-ACE  (Avva,  1993).  CFD-ACE  is  a  pressure 
based  finite-volume  code  designed  for  the  analysis  of  incompressible  and  compressible 
flowfields,  using  a  modified  SIMPLEC  algorithm.  All  solutions  used  2nd  order  accurate 
central  differencing,  stabilized  with  upwind  biasing  to  solve  the  2-d  incompressible  Navier- 
Stokes  equations.  The  computational  domain  extended  63.5  mm  (2.5  in)  upstream  of  the 
bluff  body  (x/H=-5)  and  508  mm  (20  in)  downstream  (x/H=40).  All  solutions  used  a  constant 
inflow  profile  matched  to  the  experimental  data  (see  Figure  6),  and  a  fixed  pressure  exit 
condition.  The  current  study  presents  results  from  three  solutions:  1)  a  symmetric,  steady- 
state,  2-d,  k-e  solution;  2)  a  standard  2-d,  k-e,  time-accurate  solution;  and  3)  an  RNG,  2-d,  k- 
e,  time-accurate  solution. 

The  first  solution  enforced  symmetry  by  modeling  only  half  of  the  geometry.  The 
solution  was  converged  to  steady  state  on  a  102x38-grid  (partially  shown  in  Figure  83), 
using  the  standard  k-e  model  with  wall  functions. 


Figure  83.  Grid  near  bluff  body  used  for  CFD  analyses. 
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The  Reynolds  stresses  for  the  steady-state  solution  were  calculated  from  the 


Boussinesq  equation; 


dUj 

p 

^  dxj  dx. 

where; 


(7.1) 


^,=0.09p^  (7.2) 

The  second  solution  modeled  the  entire  flowfield  on  a  102x77-grid  (mirror  of 
symmetric  grid)  and  required  a  time-accurate  solver.  Again,  the  standard  k-e  model  was  used 
with  wall  functions.  The  solution  was  first  advanced  200  iterations  in  a  steady-state  mode  to 
establish  the  flowfield,  and  also  to  allow  an  asymmetric  velocity  field  to  develop.  The 
solution  was  then  restarted  in  a  time-accurate  mode  and  advanced  4500  time-steps  at  5.0x10  * 
sec  per  time-step.  Finally  the  solution  was  advanced  -2500  time-steps  to  collect  the 
turbulence  statistics.  Care  was  taken  to  ensure  that  the  statistics  were  collected  over  a  whole 
period  of  the  vortex  shedding  cycle.  Turbulence  statistics  were  calculated  in  a  fashion  similar 
to  the  experimental  procedure,  with  each  time-step  producing  one  sample. 

The  third  solution  was  obtained  in  a  similar  manner  except  the  RNG  k-e  model  was 
employed.  The  reader  is  referred  to  Yakhot,  et  al.  (1992)  for  more  information  on  the  RNG 
theory.  One  of  the  major  differences  between  the  standard  k-e  model  and  the  RNG  k-e 
model,  is  that  many  of  the  constants  are  derived  from  the  RNG  theory,  removing  much  of  the 
ad-hoc  nature  of  the  standard  k-e  model.  The  RNG  k-e  model  was  found  by  Speziale  and 
Thangam  (1992)  to  give  good  results  for  the  backward-facing  step  flow  of  Kim,  et  al.  (1978). 
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7.3  Results  and  Discussion 


The  mean  velocities,  TJ ,  and  V,  and  three  of  the  Reynolds  stresses,  mm,  mv,  and  w 
for  all  three  solutions  are  shown  with  the  experimental  measurements  in  Figures  75,,  76,  and 
78  through  80.  Predictions  using  the  three-component  schemes  listed  above  are  also  included 
in  these  figures.  The  predicted  mean  velocities,  using  the  symmetric  steady-state  method, 
indicate  that  the  numerical  predicted  recirculation  bubble  extends  to  x/H=2,  much  further 
downstream  than  the  experimentally  determined  bubble  position.  The  numerical  prediction 
gives  a  bubble  length  of  four  half  bluff  body  heights,  similar  to  a  backward-facing  step.  This 
is  to  be  expected  since  the  bluff-body  model  resembles  a  backward-facing  step  with  a  slip 
wall  at  the  centerline  for  the  symmetric  case  method.  The  predicted  Reynolds  stresses  are 
also  comparable  to  those  measured  in  a  backward-facing  step  flows.  However,  the 
experimentally  measured  Reynolds  stresses  for  the  bluff-body  flow  were  found  to  be  ~5 
times  higher  than  those  found  in  backward-facing  step  flows.  The  underprediction  of  the 
Reynolds  stresses  by  the  steady-state  solution  leads  to  the  incorrect  bubble  length  in  the 
recirculation  zone.  For  the  purpose  of  fuel  injection  and  mixing  studies,  such  an  analysis  is 
not  sufficient. 

The  time-accurate  solutions  were  able  to  predict  the  vortex  shedding  behind  the  bluff 
body.  The  standard  k-£,  time-accurate,  solution  predicted  a  Strouhal  number,  St=0.240, 
which  compared  well  to  the  experimental  value  of  0.243.  The  predicted  time-averaged 
velocities  and  Reynolds  stresses  are  also  shown  in  Figures  75,  76,  and  78  through  80  by  the 
dot-dash  lines.  The  separation  bubble  calculated  with  the  time-accurate  solution  method  was 
considerably  shorter  than  that  predicted  using  the  steady-state  solution  method,  extending 
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only  to  x/H=1.0.  However,  the  predicted  length  of  the  bubble  was  still  longer  than  the 
measurements  indicated,  and  the  predicted  wake  spread  slightly  slower  than  the 
experimentally  determined  wake  size.  While  the  magnitudes  of  the  predicted  Reynolds 
stresses  were  higher  when  using  this  time-accurate  scheme  than  those  predicted  using  the 
steady-state  method,  they  remained  lower  than  the  experimental  results. 

To  improve  the  time-accurate  results,  the  RNG  k-e  turbulence  model  developed  by 
Yakhot,  et  al.  (1992)  was  tested.  The  Strouhal  number  predicted  by  the  time-accurate  RNG 
model  method  was  0.25,  slightly  higher  than  that  obtained  from  the  experimental  data. 
However,  the  mean  velocity  predictions  were  significantly  improved,  as  shown  in  Figures  75 
and  76  by  the  solid  lines.  The  predicted  separation  bubble  was  shorter  with  this  scheme  and 
ended  slightly  ahead  of  the  experimental  measurements  (x/H=0.7).  In  addition,  the  overall 
mean  velocity  field  was  predicted  fairly  well,  except  for  the  transverse  velocity  just 
downstream  of  the  stagnation  point.  This  deficiency  in  accurately  predicting  the  transverse 
velocities  may  be  a  result  of  3-d  effects  near  the  bluff  body. 

Comparisons  of  the  predicted  Reynolds  stresses  to  the  measurements  are  also 
improved  with  the  time-accurate  RNG  model  method  (Figure  78).  However,  the  normal 
stresses  are  overpredicted  (Figures  79,  80),  leading  to  the  shorter  separation  bubble.  The 
largest  discrepancy  occurs  with  the  transverse  normal  stress,  the  values  of  which  are 
overpredicted  by  a  factor  of  three.  Martensson,  et  al.  (1991)  obtained  similar  results  using 
LES,  and  found  that  the  transverse  normal  stresses  were  reduced  significantly  when 
simulated  with  a  3-d  analysis. 
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CHAPTERS 


CONCLUSIONS 

A  three-component  laser  Doppler  velocimeter  was  used  to  measure  the  three 
components  of  the  velocity  vector,  normal  and  shear  stresses,  and  all  nine  turbulent  triple 
products  in  the  isothermal  and  incompressible  turbulent  flowfield  downstream  of  a  2-d  bluff 
body  mounted  at  mid-stream  in  a  rectangular  test  section.  A  benchmark  quality  database  is 
established  for  future  reference  and  for  computational  fluid  dynamics  code  development  and 
validation.  The  three-component  measurements  were  validated  using  separate  one- 
component  LDV  measurements  by  showing  excellent  agreement  of  the  axial  mean  velocity 
throughout  the  tunnel.  A  mass  flow  balance  was  performed  by  numerically  integrating  the  1- 
d  data  acquired  from  top  to  bottom  of  the  tunnel.  The  agreement  for  all  measured  stations 
was  within  2%  of  the  most  upstream  station.  Negative  axial  velocities  were  found  in  the 
recirculation  zone.  Spanwise  velocity  measurements  in  the  near-field  indicated  that  the  flow 
is  three-dimensional  just  downstream  of  the  bluff  body.  The  measured  turbulent  triple 
products  are  unique  in  that  no  similar  work  has  been  reported  previously.  Validation  of  some 
commonly  used  normal  stresses,  and  turbulent  triple  product  approximations  was  made  and 
found  to  greatly  overpredict  the  experimental  values.  Newly  derived  models  for  ww,  uww, 
and  vww  in  this  flow  were  found  using  the  direct  measurements.  A  balance  of  the  turbulent 
kinetic  energy  in  this  flow  was  also  performed.  The  production,  convection,  and  diffusion  of 
turbulent  kinetic  energy  were  computed  directly  from  the  experimental  data  using  numerical 
differencing.  Turbulence  dissipation  rate  was  obtained  by  balance.  Regions  in  the  flowfield 


168 


where  turbulent  kinetic  energy  is  produced,  convected,  and  diffused  were  also  identified. 

The  measurements  also  revealed  that  the  flowfield  was  dominated  by  strong,  periodic 
vortex  shedding  in  the  wake,  yielding  a  symmetric  vortex  pair  in  the  mean  flow.  The  mean 
spanwise  velocity  component  was  measured  in  the  wake.  Significant  spanwise  flow  was 
found  at  axial  stations  less  than  x/H=1.2,  however  the  source  for  this  spanwise  component 
has  not  been  clearly  identified.  The  time-averaged  recirculation  bubble  behind  the  bluff  body 
ended  at  x/H=0.9,  significantly  shorter  than  the  recirculation  zone  found  in  a  backward¬ 
facing  step  flow.  The  dominant  vortex-shedding  frequency  was  found  to  be  454  Hz, 
corresponding  to  a  Strouhal  number  of  0.243.  The  strong  vortex  shedding  was  responsible 
for  the  presence  of  high  turbulence  in  the  near  field  with  a  high  degree  of  anisotropy. 

A  parallel  Computational  Fluid  Dynamics  (CFD)  study  was  initiated  to  gauge  the 
predictive  accuracy  of  currently  available  CFD  techniques.  Three  solutions  were  compared:  a 
2-d  steady-state  solution  using  the  standard  k-e  model,  a  standard  k-e,  time-accurate  solution, 
and  a  2-d  time-accurate  solution  using  a  Renormalized-Group  (RNG)  k-e  turbulence  model. 
The  numerical  simulations  showed  that  the  bluff-body  flowfield  cannot  be  accurately 
predicted  using  steady-state  calculations,  because  no  mechanism  exists  to  predict  the 
amplified  time-averaged  Reynolds  stresses  which  result  due  to  vortex  shedding.  The  time- 
accurate  standard  k-e  model  predicts  the  vortex-shedding  frequency  accurately,  but 
underpredicts  the  Reynolds  stresses  and  overpredicts  the  length  of  the  recirculation  zone  by 
10%.  The  RNG  k-e  model  predicts  a  higher  vortex-shedding  frequency,  but  offers  a  closer 
match  to  the  mean  velocities.  However,  the  Reynolds  stresses  tend  to  be  overpredicted.  A  3- 
d  analysis  covering  the  entire  spanwise  domain  should  improve  the  results. 
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APPENDIX  A 


LDV  Data  Reduction  Equations 

Velocity  Calculation 

The  fringe  spacing,  d,.  ,  depends  only  on  the  wavelength  of  light,  X,  and  the  full  angle 
according  to  the  following  equation; 

d  ^ 

^  2sin(<|)  /  2) 

The  velocity  is  then  obtained  from; 

U.  =  fod, 

where; 

X  =  wavelength  of  incident  light  (different  for  each  color)  in  meters, 

/„  =  frequency  detected  at  photodetector  (measured  frequency  of  particle)  in  Hz, 

=  velocity  of  flow  perpendicular  to  fringes,  in  m/sec, 

=  distance  between  fringes,  in  meters,  and 

(j)  =  angle  of  beam  intersection  (measured  prior  to  data  acquisition),  in  degrees. 

For  the  case  with  frequency  shifting,  the  following  equation  is  used; 

where; 

/,  =  shifted  frequency  (operator  preset)  =  40  MHz. 

Beam  Diameter 

The  probe  volume  diameter  equation,  based  on  l/e^  is; 


where; 


X  =  laser  beam  wavelength,  in  meter, 
f  =  lens  focal  length,  in  meter, 

D  .^  =  light  beam  diameter  entering  lens  and  defined  as  region  where  intensity  is  greater  than 
1/e^  (-0.135)  of  centerline  intensity;  1.5E-3  for  green;  1.36E-3  for  blue,  in  meter,  and 

E  =  beam  expansion  ratio,  in  this  case:  150/40=3.75. 

Probe  Volume  Length 

The  probe  volume  length,  ,  is  obtained  by; 

I  = — — - 

^  tan((t)  /  2) 


Number  of  Fringes 

The  number  of  fringes,  N ,  is  obtained  by; 


where,  d^,  and  are  obtained  as  shown  above. 
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APPENDIX  B 


Mean  Velocity  Equations  &  Coordinate  Transformation 


In  order  to  allow  for  implementation  of  nonperfectly  orthogonal  angles  as  shown  in 
illustration  B.l,  modification  of  the  equations  derived  for  orthogonal  system  was  necessary. 


Orthogonal  System 

y 


Illustration  B.l.  Departure  from  orthogonal  system. 

Where,  6^,  0^  are  the  deviation  angles  of  the  blue  and  the  green  beams  from  being  orthogonal. 
The  included  angle=9O°+0|,-6g  There  is  an  infinite  number  of  combinations  for  0^  and  0^  for 
any  configuration,  and  the  values  chosen  may  effect  8.  A  typical  use  would  be  to  input  0^=0^- 
0^  and  0^=0;  5  would  then  be  the  angle  between  the  green  beam  and  the  reference  axis. 

The  following  equations  were  used  to  perform  a  coordinate  transformation  for  an  LDV 
system  where,  the  x  &  y-coordinates  may  vary  from  being  at  90  degrees  to  each  other. 
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For  a  three-component  system,  the  following  equations  are  used: 
f/*  sinG  -t/„  sin0. 

T  7  _  ^  ^  h 

'  sin(e^-ej 
f/,,  cos6„  -U „  COS0. 

Tj _ i i _ iL 

sin(0,-0,) 

t/j  =  COS  0, 

where; 

U^,  U,,,  and  (/^  are  the  velocities  measured  by  the  tilted  system, 
t/,  ,  f/j,  and  C/j  are  the  results  of  the  transformed  velocities, 
g,  b,  and  v  are  related  to  processors  1, 2,  and  3  respectively,  and 
where; 

0^  =  the  tilted  angle  in  the  x-direction,  and 
0^=  the  tilted  angle  in  the  y-direction,  and 
0^  =  the  tilted  angle  in  the  z-direction. 


For  a  perfectly  orthogonal  system,  each  component  would  measure  directly  the  pertaining 
velocity.  But,  at  times,  even  when  great  effort  was  made  in  order  to  obtain  a  perfectly 
orthogonal  system,  it  was  found  that  this  was  very  difficult  to  accomplish.  Therefore,  these 
equations  allow  for  this  type  of  correction. 

Now,  if  the  xy-axis  is  rotated  as  was  the  case  here  (rotation  is  45  degrees),  corrections  are 
needed  that  would  allow  for  correct  measurement  of  all  three  components  (see  Dlustrations 
B.2),  where; 

6  =  the  rotation  of  the  system  in  the  xy-plane  or  about  the  z-axis  where  positive  rotation  is 
counterclockwise  (looking  from  the  left).  5=0  defines  an  LDV  system  whose  x-axis  is 
parallel  to  the  experiment  x-direction. 


181 


Positive  rotation  defined  by  Right-Hand  Rule 

Illustration  B.2.  Rotation  about  the  z-axis,  or  in  the  xy-plane. 

The  next  set  of  equations  were  used  to  perform  a  coordinate  transformation  for  rotating  a 
system  about  the  z-axis  (rotating  the  system  in  the  xy  plane)  and  rotation  of  a  system  about 
the  x-axis  (tilting  out  of  the  xy-plane,  see  Illustration  B.3). 

[/  =  [/, cos5  -(/jSinS 

V  =  C/i  cos  e  sin  5  4- [/j  cos  e  cos8  -  f/,  sin  8 
IT  =  1/,  sin  8  sin  5  +  f/j  sin  £  cos6  4- f/,  cos  8 
where; 

8  =  the  tilt  of  the  system  out  of  the  xz-plane  or  about  the  x-axis  where  positive  £  is  a  tilt 
above  the  xz-plane  and  negative  £  is  below  the  xz-plane.  8  =  0  defines  an  LDV  system  whose 
axis  is  in  the  xz-plane  (in  this  study,  £  was  set  =  0). 

U,  V,  and  W  are  the  results  of  the  transformed  velocities. 
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APPENDIX  C 


T.nv  nefinitions  and  Description  of  Components 
DeHnitions 


Laser  Doppler  velocimetry  (LDV)  is  a  technique  for  measuring  fluid  flow  velocity  by 
tagging  the  fluid  and  measuring  the  velocity  of  the  tagged  material. 

The  laser  is  a  device  for  transforming  incoherent  light  of  various  frequencies  of  vibrations 
into  a  very  narrow  intense  beam  of  coherent  light.  The  name  is  derived  from  the  initial  letters 
of  “Light  Amplification  by  Stimulated  Emission  of  Radiation.”  In  the  emission  of  ordinary 
light  the  molecules  or  the  atoms  of  the  source  emit  their  radiation  independently  of  each 
other,  and  consequently  there  is  no  definite  phase  relationship  among  the  vibrations  in  the 
resultant  beam.  The  light  is  therefore  incoherent.  The  laser,  by  means  of  an  optical  resonator, 
forces  the  atoms  of  the  material  of  the  resonator  to  radiate  in  phase.  The  emitted  radiation  is 
stimulated  by  the  excitation  of  the  atoms  to  a  higher  energy  level  by  means  supplied  to  the 
device. 

Most  lasers  consist  of  three  components,  a  resonator  structure,  a  plasma  tube,  and  a  power 
supply.  The  resonator  holds  the  two  reflective  mirrors  in  precise  alignment,  thus  forming  a 
resonant  optical  cavity.  One  of  these  mirrors  is  partially  transmissive  and  will  allow  a 
fraction  of  the  light  beam  to  escape  as  light  output  power.  The  plasma  tube  provides  optical 
gain  within  the  resonator  causing  it  to  act  as  an  oscillator.  The  power  supply  provides  d.c. 
voltage  across  the  plasma  tube  to  sustain  a  controlled  arc  discharge  through  the  fill-gas  of  the 
plasma  tube  (i.e.,  to  excite  gas  atoms  to  a  higher  energy  level  which  in  turn  radiate  to  a  lower 

state). 

Johann  Doppler  first  predicted  that  the  measured  frequency  and  wave  length  of  the  scattered 
light  depends  on  the  relative  motion  of  the  scattering  source  and  the  observer.  Now,  in  LDV 
the  particles  in  the  fluid  are  illuminated  by  the  incident  focused  laser  beams  and  become 
sources  of  scattered  light.  The  change  in  frequency  (Doppler  shift)  of  the  scattered  light  due 
to  the  motion  of  the  particles  when  observed  by  a  stationary  observer  (photomultiplier  tube) 
allows  the  particle  velocity  to  be  determined.  The  minimum  system  requirements  are: 


-  Light  source  (laser) 

-  Transmitting  optics 

-  Receiving  optics 

-  Photodetector  (converts  light  energy  to  electrical  signal) 

-  Signal  processor  (converts  frequency  to  voltage) 
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-  Data  processor  (computer  etc.):  to  calculate  moments,  store  data,  etc. 


Description  of  Components 


Laser 

The  light  source  is  usually  an  Argon-ion  laser  which  emits  strongly  at  476.5, 488,  514.5  nm. 
For  most  four-beam  two-color  systems,  the  488  (blue)  and  514.5  (green)  nm  lines  are  used. 

Transmitting  Optics 

1.  Beam  collimator  -  is  usually  mounted  directly  on  the  laser  output  port.  This 
component  ensures  minimum  beam  divergence  and  locates  the  beam  waist  at  the  focal  point 
of  the  final  focusing  lens. 

2.  Color  separator  -  This  component  separates  the  output  laser  beam  into  individual 
colors  by  means  of  a  dispersion  prism  assembly. 

3.  Plane  mirrors  -  These  serve  to  align  the  individual  beams  with  the  previously  defined 
optical  axis. 

4.  Polarization  rotators  -  These  components  ensure  correct  orientation  of  the  incident 
beam  polarity. 

5.  Beam  splitters  -  These  components  split  the  incident  beam  into  two  beams  of  equal 
intensity. 

6.  Bragg  cells  -  These  components  are  acousto-optic  devices  which  are  used  to 
frequency  shift  one  of  the  beams  from  the  beam  splitter.  The  other  beam  is  passed  through  an 
optical  block  to  maintain  coherence. 

7.  Beam  steering  modules  -  These  are  small  prisms  which  are  used  to  steer  the  beams 
for  proper  alignment  and  ensuring  the  beams  cross  at  the  same  point. 

8.  Beam  spacers  -  These  components  serve  to  reduce  the  spacing  between  the  parallel 
beams  before  entering  the  beam  expander. 

9.  Beam  expander  -  This  component  serves  to  increase  the  beam  diameter  and  spacing, 
thereby  resulting  in  a  finer  probe  volume. 

10.  Focusing  lens  -  This  component  focuses  the  beams  from  the  beam  expander  to  form 
the  probe  (measurement)  volumes. 
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Receiving  Optics 


1.  In  the  back-scatter  mode  the  focusing  lens  acts  as  the  collector  of  the  scattered  light 
from  the  probe  volume. 

2.  Field  stop  system  -  This  is  used  to  optimize  the  focal  depth  of  the  collecting  optics, 
i.e.  it  serves  as  an  optical  filter  for  the  elimination  of  extraneous  background  light. 

3.  Color  separator  -  In  a  two  color  system,  this  component  separates  and  directs  the 

scattered  light  by  means  of  a  dichroic  mirror. 

4.  Receiving  modules  -  These  components  consist  of  a  lens  and  mirror  assembly  which 
focuses  the  scattered  light  on  to  the  photomultiplier  tube. 


APPENDIX  D 


Data  Reduction  Notes 


Mean: 


RMS: 


U  =  -tu,  or 

n  '■=' 


f/™  = 


SIGMA:  C/„  = 


p; 

t'n 

1=1 

? - 

_ 

i  ' 

n 

1  ) 

or 


or 


+ul+u;+....ul 

n 


uL -u‘ 


where  n  is  the  number  of  samples  in  the  data  set. 

Turbulence  intensity(%):  100  Ua  /U^, 

Second  Moments:  «v  =  —  {{u.  -  U )(v,  -  V )) 

=-t{u,v,-u,v-uv,+vv) 

n  '■=!  J  y  M  n  J  \  n  J 
ri  ^  _ _ _ 

-tuy.  -(vf/)-(t/v)+(c/y) 


=  (-tuy^-{vu)=uv  -uv 

yn  -=i  j 


i  /=1  W  y  V  i=i  W  A  1=1 


i=l 


Similarly, 

Triple  Products: 


''  U.U.'^  ^  ^  t7  V 


uu  = 


t- 


/=i  n 


\^>=i  n  n 


uvw 


=^(t{u,-uXv,-v)(w,-w)) 

=  -^(t(f/,y, -u.v-uv,.+uv){w,  -ir)) 
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vw.  -u.vw-  uvw. + UVW.  -UV.W+V  VW+UVW-  uvw 


ill  I  I  ft 


=  f- i  1-fv 


w- tf/y  l+[ywi^  l-(f/vi^) 

« -•=!  '  '  J  {  '■='  «  y  V  '='  n  J 


/=!  n 


u-tvw..  =  ut- 


Uvt^  =  UVW;  W-luy,=wt- 


Vwt—  =  VWU; 

i=i  n 


uwl—=uwv 


j^u}W_yf^^l^_jjf^YK+uyw_^^±^+vwu+uwv-uvw 


/=!  n 


=1  n  i=i  n 


=1  n  '='  n 


Similarly: 


uuv  =  ^{i(u.-U)  (v,-vfj 

^^{t{u,U,-2U,U  +  UU){v,  -V)) 


=-  I(t/,c/,y, -u^uy-2u,uv.+2u.uv+uuv^ -uuv, 

^  _2^j^^+2(t7Z7\7) 

/=!  n  '=•  n  '=*  w 


Similarly; 


UUW.  UJJ 


;=l  n 


-wl^-2ut^+2{UUW) 


■=±Wl^U±YK-2vt^+2{UVV) 


i=\  n  '='  u  '='  n 


^=t^^^-wt^-2vt^+2{VVW) 

i=i  n  i=i  n  <=i  n 


188 


1=1  n  '=1  n  /=!  n 

f=i  w  /=i  n  /=!  n 

=  i  ^  -  3Z7i  + 2t7' 

i=i  n  1=1  n 

__  n  v^  — "  yy.  —1 

wv  =  -  3yl-^ + 2y 

i=i  n  1=1  n 

-  — aW.  —3 

WW  =  I-^ — 3Wl-^+2W 
/=i  n  1=1  n 

Skewness  Coefficient;  3"*  Moment/2aJ 

4"’  Moment  (U) . t^-4ut’^+6U'tuf -3U' 

t=i  n  '=!  n  »’=i 

4"'  Moment  (V) . t- — 4VX^+  ev^tv^  -  3^ 

/=!  n  M  n  i=i 

«  W'^  _ n  _ .  n  _ . 

4"’  Moment  (W) . Z-^-4\yI-^+  6iy  IW.'  -3iy 

,=i  n  1=1  n  t=i 

Flatness  Coefficient  (Kurtosis);  4*^  Moment/a 

Note: 

The  3"'  Moments  are  also  called  Skewness;  and  the  4"’  Moments  are  called  Flatness. 


% 


r 
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APPENDIX  E 


Software  File  Output  Format 


X(in)  Y(in)  Z(in) 

Uref(in/s) 

U/Uref 

V/Uref 

W/Uref 

Usig/Uref 

Vsig/Uref 

Wsig/Uref 

uulU^ref 

wlU^ref 

^lU^ref 

uv  1  U^ref 

uw/U^ref 

ml  U"- ref 

mu  /  U^ref 

uuv  /  U^ref 

uww!  U^ref 

vuu/  ref 

'^/U^ref 

vvwl  U^ref 

wuu  /  U^ref 

vww/U^ref 

WWW  !  ref 

r 
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APPENDIX  F 


Turbulent  Kinetic  Energy  Equation 

Starting  with  the  Navier  Stokes  Equation  in  tensor  notations,  without  the  time-dependent 
term; 


pUj 


Substitute; 


U.=U.  +  u.  &  P  =  P+p 
to  get; 


p(f/  i  +  Uj) 


Bx. 


B(P+P)  d\Ui  +  u.) 
+  ji 


Bx 
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Now  multiply  each  term  by  u.  and  take  the  mean  of  this  scalar  equation,  realizing  that  the 
product  of  a  mean  and  a  fluctuating  quantity  has  zero  mean;  and  the  mean  of  a  fluctuating 
quantity  is  zero  as  well. 
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Note  that; 
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Therefore;  the  equation  reduces  to; 
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Also,  the  following  modifications  may  be  made; 
—  Bu.  —  B  uu. 
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since;  7—^  =  — ^  =  0 ,  from  the  continuity  equation. 
dx,  dx. 


After  these  substitutions  are  made,  and  letting  k  =  -(m,m,),  we  get; 
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or  in  another  form; 
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Where; 
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pU  - —  =  Convection  terms 
^dx. 


1  3 


3(m,p) 


- ( pu  u.u  )-— nr —  =  Diffusion  terms 

2  3x.^‘^^''  dx. 


- ^  . 

-M  M .  -T —  =  Production  terms 

'  ^  dx. 


d^u 


+|IM,  =  Dissipation  terms 

OX 


(F.8) 

(F.9) 

(F.IO) 


(F.il) 

(F.12) 

(F.13) 

(F.14) 

(F.15) 

(F.16) 


192 


'I 


r 


Of  this  equation,  the  following  terms  were  estimated  using  experimental  measurements; 

1 .  All  convection  terms 

2.  All  production  terms 

3.  All  gradients  of  triple  products 


Dissipation  was  calculated  by  balance  assuming  the  terms  involving  pressure  are  negligible. 
However,  dissipation  (per  unit  mass)  is  defined  as  (Wilcox,  C.  1993); 


E  =-v 
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where,  this  term  is  obtained  from; 
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a  dk 

- — (v  r — )  =  Molecular  Diffusion 
ox.  Bx. 


If  molecular  diffusion  is  negligible,  than  e  is  equivalent  to  the  dissipation  terms  shown 
above. 


The  turbulent  kinetic  energy  equation  is  similar  to  the  Reynolds-stress  equation  (shown  in 
the  next  page)  with  the  exception  of  the  missing  Pressure-Strain  interaction  term.  The  reason 
for  this  may  be  deduced  from  Equation  (9)  above; 
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which  shows  that; 
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to  satisfy  mass  continuity. 

Therefore,  the  Pressure-Strain  interaction  terms  make  no  contribution  to  the  overall  level  of 
turbulence  activity. 
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Reynolds  Stress  Transport  Eanatton 


d  u.Uj 

dx, 
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Of  this  equation,  the  following  terms  may  be  experimentally  estimated  using  experimental 
measurements: 


1.  All  convection  of  Reynolds  stresses  terms 

2.  All  production  of  Reynolds  stresses  terms 

3.  All  gradients  of  triple  products 

The  pressure-diffusion  terms  are  very  difficult  if  not  impossible  to  measure  at  this  point  and 
are  usually  assumed  negligible.  This  assumption  is  reasonable  for  fully  turbulent  flows  at 
locations  away  from  the  wall,  such  as  is  the  case  here.  The  pressure-strain  interaction  terms 
are  also  referred  to  as  the  ‘redistribution’  terms  because  they  do  not  contribute  in  the 
generation  or  the  destruction  of  turbulence,  but  they  simply  smooth  the  distribution  of  the 
stresses  throughout  the  flow.  This  term  vanishes  for  constant  density  flows  because  of  the 
mass  continuity  constraint  (Hogg,  S.,  1 988): 


=  0 


The  viscous  dissipation  terms  contain  the  dissipation  of  the  Reynolds  stresses  due  to  small 
scale  molecular  motion  and  interaction. 
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Uncertainty  Calculation 


Bernoulli’s  equation:  AP  =  — 


Density  ( p )  = 


RT 


Velocity  (C/)  = 


2AP 


K  RT, 


U„,  (max)= 


i 


2AP^ 


rt: 


U  .  (min)= 


2AP- 


|pr; 


Humiditv  Correction  Error  to 

Since  the  one-component  data  were  taken  in  winter,  while  the  3-component  data  in  the 
summer,  the  reference  velocity  error  due  to  the  changes  in  the  relative  humidity  would  be 
computed  as  follows; 


R=RI  M.W. 


where;  M.W.  =  M.W.  (air)  +  M.W.  (air)  x  05 
M.W.  =  molecular  weight 
03  =  humidity  ratio  =  0.622  x  Pv/Pa 
Pv  =  O  Pg 

Pg  =  saturation  pressure  at  P 
<I>  =  relative  humidity 
Pa  =  P  -  Pv 


P  =  atmospheric  pressure 
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AC7  =  ±1.96CT„/Viv 
AV  =  ±l96aJyfN 
AW  =  ±\96^j4n 


Uncertainty  in  the  Means 


where;  the  constant  1.96  corresponds  to  a  confidence  interval  of  95%,  and  a„,  and  are 
the  estimators  of  the  largest  standard  of  deviation  values  of  the  U,  V,  and  W-components. 

N  is  the  number  of  samples. 

To  obtain  the  percent  error,  AU  is  then  divided  by  ,  the  upstream  reference  velocity,  the 
result  is  subtracted  from  1 ;  and  the  outcome  is  multiplied  by  100. 

(Note:  for  a  99%  confidence  interval,  the  constant  to  be  used  is  2.576  (instead  of  1.96). 


APPENDIX  H 


Uncertainties  in  TKE  (k)  and  TKE  Terms 


The  turbulent  kinetic  energy  “  A:  ”  is  the  sum  of  three  terms  and  is  written  as  such; 
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k  =  —(uu  +  vv  +  ww) 


(H.1) 


In  order  to  obtain  the  statistical  uncertainty  for  “  k  ”,  the  “Sum”  rule  needs  to  be  applied. 

First,  recall  that  the  basic  measure  of  uncertainty  is  the  variance  (the  actual  values  to  be  used 
are  those  of  the  uncertainty  errors  obtained  with  the  “jackknife”  method  described  by 
Benedict  &  Gould  (1995)  and  are  listed  after  each  parameter  in  the  computer  output).  There 
are  basic  rules  for  linear  combinations  of  random  variables.  For  instance: 

If  a  and  b  are  constants,  and  X  and  Y  are  random  variables,  then; 
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If  X  and  Y  are  independent,  than  the  last  term  drops  out; 


(H.4) 


The  variance  for  a  difference  (when  X  and  Y  are  independent)  is  the  same  as  the  sum; 


^2  _  „2_  2  ,  ,  2_  2 
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To  get  the  uncertainty  in  any  of  these  terms,  the  square  root  needs  to  be  calculated. 
As  an  example  of  a  “Sum”  term,  the  uncertainty  for  “k  ”  is  evaluated  as  such; 

1  .  1 


1  ,2  1 


+T<Jw  =T(fyi 


(H.6) 


and; 


(H.7) 


Therefore,  the  uncertainty  in  “A”,  is  the  square  root  of  the  sum  of  the  square  of  the 
associated  uncertainties  of  each  term. 
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Note:  there  is  also  another  method  to  calculate  the  uncertainty  in  a  “Sum”.  The  error 
equation  is; 


=  (H.8) 

where,  0;  are  the  partial  derivatives,  and  Gj  are  the  variances.  Using  this  equation  to  find  the 
uncertainty  in  “  A:  ”  gives: 


ft  .  ft  =  1  .  ft  =  1 

'  duu  2  ’  ^  aw  2  ’  ^  dww  2 


and. 


(H.9) 


(H.10) 


<7*  =  + \^L  =  l^l^uu+<  +<^L  (sameasH.7)  (H.ll) 

Now,  for  the  uncertainty  in  a  “Product”  term,  equation  H.8  still  applies.  For  instance,  the 
uncertainty  in  the  “uu/  k  ”  term  is  calculated  as  such; 

let:  P  =  ^  (H.12) 

k 


then: 
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1 


For  the  “TKE”  terms,  the  same  technique  is  used.  However,  now  one  has  to  deal  with  “A:  ”  in 
addition.  For  instance,  in  the  “Convection”  of  TKE  term,  Udk  /  dx ,  the  following  is  u.sed; 
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The  first  term  on  the  left  of  (H.16)  is  obtained  from  the  experimental  data,  using  central 
differencing.  The  second  term  is  the  uncertainty  in  the  mean  axial  velocity,  and  is  available. 
The  third  term  is  the  mean  axial  velocity.  The  last  term  in  (H.16)  is  dealt  with  as  follows; 


(H.17) 


and  the  terms  in  the  numerator  are  calculated  according  to; 


(H.18) 


and. 


(H.19) 


Similar  calculations  need  to  be  done  for  Vdk  /  dy ,  and  for  all  other  TKE  terms. 
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